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i— i. Abstract 

Q-f We study here dissipative properties of systems composed of two components one of 

which is highly lossy and the other is lossless. A principal result of our studies is that 
all the eigenmodes of such a system split into two distinct classes characterized as high- 
loss and low-loss. Interestingly, this splitting is more pronounced the higher the loss of 
the lossy component. In addition, the real frequencies of the high- loss eigenmodes can 
become very small and even can vanish entirely, which is the case of overdamping. 

> 

1 Introduction 
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We introduce a general framework to study dissipative properties of two component systems 
composed of a high-loss and lossless components. This framework covers conceptually any 
dissipative physical system governed by a linear evolution equation. Such systems include in 
particular damped mechanical systems, electric networks or any general Lagrangian system 
with losses accounted by the Rayleigh dissipative function, [Pars} Sec. 10.11, 10.12], |Gant| 
Sec. 8, 9, 46]. 

An important motivation and guiding examples for our studies come from two-component 
dielectric media composed of a high-loss and lossless components. Any dielectric medium 
always absorbs a certain amount of electromagnetic energy, a phenomenon which is often 
referred to as loss. When it comes to the design of devices utilizing dielectric properties very 
often a component which carries a useful property, for instance, magnetism is also lossy. So 
the question stands: is it possible to design a composite material/system which would have a 
desired property comparable with a naturally occurring bulk substance but with significantly 
reduced losses. In a search of such a low-loss composite it is appropriate to assume that 
the lossy component, for instance magnetic, constitutes the significant fraction which carries 
the desired property. But then it is far from clear whether a significant loss reduction is 
achievable at all. It is quite remarkable that the answer to the above question is affirmative, 
and an example of a simple layered structure having magnetic properties comparable with a 
natural bulk material but with 100 times lesser losses in wide frequency range is constructed 



in FigVit8 . The primary goal of this paper is to find out and explain when and how a two 
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component system involving a high-loss component can have low loss for a wide frequency 
range. The next question of how this low loss performance for wide frequency range is 
combined with a useful property associated with the lossy component is left for another 
forthcoming paper. 

A principal result of our studies here is that a two component system involving a high-loss 
component can be significantly low loss in a wide frequency range provided, to some surprise, 
that the lossy component is sufficiently lossy. An explanation for this phenomenon is that if 
the lossy part of the system has losses exceeding a critical value it goes into essentially an 
overdamping regime, that is a regime with no oscillatory motion. In fact, we think that for 
Lagrangian systems there will be exactly an overdamping regime but these studies will also 
be conducted in already mentioned forthcoming paper. 

The rest of the paper is organized as follows. The model setup and discussion of main 
results and their physical significance is provided in Section |2j In the following Section [3] we 
apply the developed general approach to an electric circuit showing all key features of the 
method. Sections H] and is devoted to a precise formulation of all significant results in 
the form of Theorems, Propositions and so on. Finally, in the last Section E] we provide the 
proofs of these results. 

2 Model setup and discussion of main results 



A framework to our studies of energy dissipation is as in |FigSchl| and |FigSch2| . This 
framework covers many physical systems including dielectric, elastic and acoustic media. Our 
primary subject is a linear system (medium) whose state is described by a time dependent 
generalized velocity v (t) taking values in a Hilbert space H with scalar product (•,•). The 
evolution of v is governed by a linear equation incorporating retarded friction 

POO 

md t v (t) = -iAv (t) - / a(r)v(t - r) dr + f (t) , (1) 

Jo 

where m > is a positive mass operator in H, A is a self-adjoint operator in H, f (t) is a 
time dependent external generalized force, and a(t), t > 0, is an operator valued function 
which we call the operator valued friction retardation function |KubTod2~| Section 1.6], or just 
the friction function. The names "generalized velocity" and "generalized force" are justified 
when we interpret the real part of the scalar product Re {(v (t) , / (t))} as the work W done 
by / (t) per unit of time at instant t, that is 

POO 

W= / Re{(v(t)J(t))} dt. (2) 



The internal energy of the system is naturally defined as 

internal energy U = — (v(t), mv(t)) , (3) 

and it readily follows from (TT]) that it satisfies the following energy balance equation 
dU 



df Re{(v(t),f(t))}-Re^(v(t), I a(r)v (t - r) dr ) \ . (4) 
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The second term of the right hand side of (jl]) is interpreted as the instantaneous rate of "work 
done by the system on itself," or more properly the negative rate of energy dissipation due 
to friction. 

If we rescale the variables according to the formulas 

mv, tt = -^=A-^=, a = -^a-U, / = — U/, (5) 
m Jm \ m ^/m Jm 



then the equation (DO) reduces to the special form in which m in the new variables is the 
identity operator, i.e. 

d t v (t) = -iQv (t) - / a (r) v (t - r) dr + f (t) , (6) 
Jo 

and in view of ([3]) the internal energy U associated with the state v turns into the scalar 
product, that is 

internal energy U = — (v(t),v(t)) . (7) 

The system evolution equation of the special form ([6]) has two important characteristic prop- 
erties: (i) the operator Q = — j=yl^= can be interpreted as the system frequency operator in 
H; (ii) the system internal energy is simply the scalar product (JTJ). We refer to this important 
form ()6]) as the system canonical evolution equation. 

For the sake of simplicity we assume the friction function to be "instantaneous," that is 
a (t) = (3B5 (t) with B self-adjoint and /3 > is a dimensionless loss parameter which scales 
the intensity of dissipation. Of course, such an idealized friction function is a simplification 
that we take to avoid significant difficulties associated with more realistic friction functions 
as in 



FigSchl , FigSch2 



Now, assuming that the rescaling flHj) was applied, the canonical evolution equation (J6J) 
takes the form 

3 t v it) = -L4 (/3) v(t) + f (t) , where A (ft) = Q - i/3B, f3>0. (8) 

Importantly, for the system operator Q — i/3B it is assumed the operator B satisfies the power 
dissipation condition 

B>0. (9) 
The general energy balance equation (jl]) takes now a simpler form 

^L = Re{(v(t)J(t))}-W dis . (10) 

where 

U = system energy = - (v(t),v(t)) , (11) 
Wdis = system dissipated power = (3 (v (t) , Bv (t) ) . 

Most of the time it is assumed that the system governed by (jHJ) is at rest for all negative 
times, i.e. 

v(t) = 0, f(t) = 0, t<0. (12) 
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To simplify further technical aspects of our studies and to avoid the nontrivial subtleties 
involved in considering unbounded operators in infinite-dimensional Hilbert spaces, we as- 
sume the Hilbert space H to be of a finite dimension N. Keeping in mind our motivations 
we associate the operator B with the lossy component of a composite and express the signif- 
icance of the lossy component in terms of the rank N B of the operator B. Continuing this 
line of thought we introduce a space H B , the range of the operator B, and the corresponding 
orthogonal projection Pb on it, that is 

H B = {Bu :ueH}, N B = dimH B . (13) 

In what follows we refer to H B as a subspace of degrees of freedom susceptible to losses or just 
loss subspace. We also refer to the orthogonal complement of = H H B as the subspace 
of lossless degrees of freedom or no-loss subspace. The number N B plays an important role 
in the Livsic theory of open systems where N B is called "the index of non-Hermiticity" of 
the operator Q — i/3B, |Livl pp. 24, 27-28]. The definition of P B readily implies the following 
identity 

B = P B BP B . (14) 
We suppose the dimension N B to satisfy the following loss fraction condition 

N R 

0<6 B = ^<1, (15) 

which signifies in a rough form that only a fraction 5 B < 1 of the degrees of freedom are 
susceptible to lossy behavior. As we will see later when the loss parameter |3 > 1, only a 
fraction of the system eigenmodes are associated with high losses and this fraction is exactly 
5 B . For this reason we may refer to 5 B as the fraction of high-loss eigenmodes. 

It turns out that the system dissipative behavior is qualitatively different when the loss 
parameter (3 is small or large. It seems that common intuition about losses is associated with 
the small values of 0. The spectral analysis of the system operator Q — i/3B for small (3 can 
be handled by the standard perturbation theory, |Bau85] . The results of our analysis are 
contained in Theorem [15] and may be summarized as follows: Let Uj, 1 < j < N denote the 
all eigenvalues of the operator Q repeated according to their multiplicities then there exists 
a corresponding orthonormal basis of eigenvectors Uj, 1 < j < N such that if < <C 1 
then the operator Q — ij3B is diagonalizable with a complete set of Q ((3) eigenvalues and 
eigenvectors Vj ((3) having the expansions 

Q ((3) = u 3 - ifi (u 3 , Buj) +0{(3), to, = { Uj , Q Uj ) , v 3 {/3) = Uj + O (/?) , 1 < j < N, < 1. 

(16) 

The effect of small losses described by the above formula is well known, of course, see, for 
instance, |Gant| Sec. 46]. 

The perturbation analysis of the system operator Q — i/3B for large values of the loss 
parameter (3 ^> 1 requires more efforts and its results are quite surprising. It shows, in 
particular, that all the eigenmodes split into two distinct classes according to their dissipative 
behavior : high-loss and low-loss modes. We refer to such a splitting as modal dichotomy. 

In view of the above discussion we decompose the Hilbert space H into the direct sum of 
invariant subspaces of the operator B > 0, that is, 

H = H B @H^ (17) 
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where Hb = ranP is the loss subspace of dimension N B with orthogonal projection P B and 
its orthogonal complement, Hg = ker B, is the no- loss subspace of dimension N — Nb with 
orthogonal projection Pg. Then the operators Q and B, with respect to this direct sum, are 
2x2 block operator matrices 



" n 2 


e " 


, B = 


' B 2 


" 


e* 











(18) 



where Vt 2 : = Pb^Pb\h b '■ Hb H B and B 2 := PbBP b \ Hb : H B — > H B are restrictions of the 
operators f2 and B respectively to loss subspace H B whereas Vt\ := PaQP B \ ± : Hr — > Hi; 

is the restriction of Q to complementary subspace Hg. Also, 6 : Hg — > P B is the operator 
6 := P B nP^L ± whose adjoint is given by 6* = P^QP B \„ : H B -»■ P^. The block 

I P*£J ' "S 

representation ffTS]) plays an important role in our analysis involving the perturbation theory 
as well as the Schur complement concept described in Appendix [7J 



2.1 Modal dichotomy for the high-loss regime 

Notice first that in view of (|S]) the operator (— if3)~ 1 A (f3) = P + i/3 _1 fi is analytic in in 
a vicinity of (3 = oo. Let then £ ((3) be an analytic in (3~ x eigenvalue of A in the same 
vicinity, with the possible exception of a pole at — oo. Notice that if use the substitution 
e = (— \(3)~ l the operator ^(ie -1 ) = P + eQ is a self-adjoint for real e and consequently 
the eigenvalue A (e) = e( (i£ _1 ) of the operator B + eQ must be an analytic function of e 
and real- valued for real e. Hence it satisfies the identity A (e) = A (e) where e is the complex 
conjugate to e. The later in view of the identity ( (j3) = (— i(3) A ((— i/3) -1 ) readily implies 
the following identities for the eigenvalue ( (/3) for real (3 in a vicinity of (3 = oo 

CC9) = CH3), or ReC(-/3)=ReC(/3), Im( (-/?) = - Im( (/?) • (19) 

Consequently, Re£ (/3) and Im£ {(3) are respectively an even and an odd function for real f3 
in a vicinity of f3 = oo implying that their Laurent series in (3" 1 have respectively only even 
and odd powers. 

The perturbation analysis for (3 ^> 1 of the operator A ((3) = Q — if3B described in Section 
14. 1| introduces an orthonormal basis {vjj}^ =l diagonalizing the operators Qi and B 2 from the 
block form (JIB"]) , that is 

B 2 Wj = Qwj for 1 < j < N B ; MiWj = pjWj for N B + 1 < j < N, (20) 

where 

Q = {w j} B 2 Wj) = (w v Bw 3 ) for 1 < j < N B ; (21) 
Pj = QiWj) = (ibj, Qwj) for N B + 1 < j < N. 

The summary of the perturbation analysis for the high- loss regime /3 ^> 1, as described in 
Theorem is as follows. The system operator A(j3) is diagonalizable and there exists a 
complete set of eigenvalues Q {(3) and eigenvectors Wj (j3) satisfying 

A(P)w j (P) = Q(P)w j (P), l<j<N, /3»1 (22) 
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which split into two distinct classes 

high-loss: Q (0) , Wj (0), 1 < 3 < Nb; (23) 
low-loss: Q (0) , (ft), N B + l<j< N, 

with the following properties. 

In the high-loss case the eigenvalues have poles at ft = oo whereas their eigenvectors 
are analytic at ft = oo, having the asymptotic expansions 

Cj{j0) = -Kj0 + Pj + O(0- X ) 1 (j>0, Pj eR, w j (ft)=w j + 0(ft- 1 ), l<j<N B . 

(24) 

The vectors Wj, 1 < j < N B form an orthonormal basis of the loss subspace Hb and 

Bwj = CjWj, pj = (wj, flwj) , for 1 < j < N B . (25) 

In the low-loss case the eigenvalues and eigenvectors are analytic at ft = oo, having 
the asymptotic expansions 

Q (ft) = p. - idjft' 1 + O (ft- 2 ) , pj E R, dj > 0, (26) 
Wj (ft) = Wj + w^ft- 1 + O (ft~ 2 ) , N B + l<j<N. 

The vectors Wj, Nb + 1 < j < N form an orthonormal basis of the no- loss subspace and 

Bwj = 0, pj = (wj, ttwj) , dj = Bw { r 1] ^ for N B + l<]<N. (27) 

The expansions (1241) and (T2"6"|) together with (1191) readily imply the following asymptotic 
formulas for the real and imaginary parts of the complex eigenvalues Cj(ft) 101 ft ^ 1 

high-loss: ReQ(ft) = p,+0(ft- 2 ), ImQ (ft) = -Q0 + O (ft- 1 ) , 1 < j < N B ; (28) 

low-loss: Re Q(ft) = Pj + O (ft- 2 ) , Im Q (ft) = -djft' 1 + O (ft~ 3 ) , N B + 1 < j < N. 

(29) 

Observe that the expansions fl28|) and fl29l) readily yield 

lim ImQ (ft) = -oo for 1 < j < N B ] lim ImQ (ft) = for N B + 1 < j < N, (30) 

justifying the names high-loss and low-loss. Notice also that the relations (T24|) — (T2TT) imply 
that the high-loss eigenmodes projection on the no-loss subspace is of order ft~ l in 
contrast to the low-loss eigenmodes for which the projection on the loss subspace Hb is of 
order ft~ l . In other words, for ft ^> 1 the high-loss eigenmodes are essentially confined to the 
loss subspace Hb whereas the low-loss modes are essentially expelled from it. 

2.2 Losses and the quality factor associated with the eigenmodes 

Here we consider the energy dissipation associated with high-loss and low-loss eigenmodes. 
The power dissipation is commonly quantified by the so called quality factor Q that can 
naturally be introduced in a few not entirely equivalent ways, [Paint pp. 47, 70, 71]. The most 
common way to define the quality factor is based on relative rate of the energy dissipation 



6 



per cycle when the system is in a state of damped harmonic oscillations v (t) with a given 
frequency u>, namely, 



^ ^ energy stored in system [/ (t>(t), t>(t) ) 

energy lost per cycle Wdis 2/3 -Bv (t)) ' 



where we used for the system energy U and the dissipated power Wd; s their expressions (TTTj) . 
Notice also that in the above formula we use the absolute value \u\ of the frequency u; since 
in our settings the frequency u can be negative. The state of damped harmonic oscillations 
v(t) is defined by an eigenvector w of the system operator A = Q — i/3B with eigenvalue £, 
and it evolves as v(t) = we" 1 ^ with the frequency u = Ke( and the damping factor — lm(. 
Its system energy U, dissipated power Wdis, and quality factor Q satisfy 



where 



U = U[w}e 2I ^ t , W dis = W dis He 2ImC ', Q = Q[w], (32) 

U[w] = ^ {w, w) , ^is H = -2 Im (U [w] , (33) 

(w,Qw) (w,f3Bw) 

ReC=^ r , ImC = 7 r-, (34) 

[w,w) [w,w) 

0M = -^- (35) 

For an eigenvector w we refer to the terms U [w], Wdis [w], and Q [w] as its energy, power 
of energy dissipation, and quality factor, respectively. Observe, that eigenvectors with the 
same eigenvalue have equal quality factors. Notice also that an eigenvector w with eigenvalue 
C has power of energy dissipation Wdi S [w] equal to the product — (w,w)lm( and quality 
factor Q[w] equal the ratio J^j^ ■ 

Consider now the high-loss regime ^> 1. Let Q ((3), 1 < j < N denote the high-loss 
and low- loss eigenvalues of the system operator A (/3) which have the expansions ( 128|) . ( |29|) . 
Then for any eigenvectors Wj (/3), 1 < j < N with these eigenvalues, respectively, which are 
normalized in the sense 

( Wj (f3) , w 3 (/?)) = 1 + (/3- 1 ) , 1 < j < N (36) 

as (3 — » oo, the following asymptotic formulas holds as /3 — > oo for the energy and the power 
of energy dissipation of these modes 

U [wj (/?)] = \ + (/T 1 ) , l<j<iV; (37) 

high-loss: W dis K (/3)] =4/3 + 0(1), 1<J<N B ; (38) 
low-loss: Wdi S [wj (/?)] = djfi- 1 + O (/T 2 ) , A^ B + 1 < j < N. (39) 

We see clearly now the modal dichotomy, i.e. eigenmode splitting according to their dissipa- 
tive properties: high-loss modes Wj (f3), 1 < j < N B and low- loss modes Wj (/3), N B + 1 < 
j < N. Indeed, these asymptotic formulas (!38l) . (139|1 imply 

high-loss modes: lim Wa s [wj (/?)] = oo; low-loss modes: lim Wdis [ w j {P)} — 0. (40) 
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The quality factor Q [wj (/?)] for each high-loss eigenmode has a series expansion con- 
taining only odd powers of (3~ l with the asymptotic formula as ft — » oo 

Q K (P)] = ^P' 1 + o (/r 3 ) , i<j<iv B . (4i) 

The quality factor Q [wj (/?)] for each low-loss eigenvectors has a series expansion containing 
only odd powers of (3~ x as well provided IvaQ ((3) ^ for (3 3> 1. Moreover, it satisfies the 
following asymptotic formula as (3 — )■ oo 

Q K- 09)] = + O (r 1 ) , iV fl + l<j<iV, (42) 

provided dj ^ 0. In fact, it is true under rather general conditions that dj > 0, for j = 
Nb + 1, . . • , N (see (1261) and Remark [9]). These asymptotic formulas (]4T1) . ( l4"2l) readily imply 
that 

high-loss modes: lim Q [wj (/?)] = 0; low-loss modes: lim Q [wj (/?)] = i °° *r ^ J ^ n' 

/3— s>oo /3— >oo I U II Pj — U. 

(43) 

Observe, that the relations and (f^^j clearly indicate that for the low-loss modes the larger 
values of (3 imply lesser losses and the possibility of a higher quality factor! In particular, 
the more lossy is the lossy component the less lossy are the low-loss modes. This somewhat 
paradoxical conclusion can be explained by the fact that the low-loss eigenmodes are being 
expelled from the loss subspace Hb in the sense that their projection onto this subspace 
satisfies asymptotically PsWj {(3) = 0(/3 _1 ) as (3 — > oo. 

2.3 Losses for external harmonic forces 

Let us subject now our system to a harmonic external force / (u) e~ luJt of a frequency u that 
will set the system into a stationary oscillatory motion of the form v (u) e~ lut of the same 
frequency u and amplitude v (u) depending on the energy dissipation. Or more generally we 
can subject the system to an external force / (t) and observe its response v (t) governed by 
the evolution equation ([S]). The solution to this problem in view of the rest condition f[T2"j) 
can be obtained with the help of the Fourier-Laplace transform 

POO 

e l ^v{t)dt, Im£ > 0, (44) 

Jo 

applied to the evolution equation ([HD resulting in the following equation 

££(£) = [Q-i/35]£(0+i/(£), f = u + iri, r/ = Im£ > 0. (45) 
For Im£ > in view of B > the operator £1 — (Q — i(3B) is invertible, and hence 

*(0 = a(0/(0> e = ^ + ir/, r/ = lme>0, (46) 
2l(e) = i[e/-(^-i/?S)]- 1 . (47) 

For a harmonic force / (t) = fe~ luJt the corresponding harmonic solution is v (t) = ve~ lwt 
where 

v = f&(u>)f, = i[ul- A(/3)] _1 , A(f3) = O - i0B, > 0. (48) 
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The operator 21 (cu) is called the admittance operator. 

For the stationary regime associated with a harmonic external force / (t) = fe~ luJt the 
quality factor can be naturally defined by a formula analogous to (13T|) . namely 



Q = Qf,u = 27T 



energy stored in system 
energy lost per cycle 



(49) 



where the energy lost refers specifically to the energy loss due to friction in the system. By 
the expressions (fTTj) . the quality factor turns into 



Q 



\CU\ 



U 



\cu\ 



/3(v,Bv) : 



(50) 



where according to (|T0|) 



[/ = - (v, v) is the stored energy, 
Wdis = /3 (v, Bv) is the power of dissipated energy 



(51) 



In many cases of interest the external force / is outside the loss subspace corresponding 
to a situation when the driving forces/sources are located outside the lossy component of 
the system. This important factor is described by the projection on the no-loss space Hg = 
HqHb, that is by Pgf. We may expect the effect of losses to depend significantly on whether 
Pgf = or Pgf 7^ 0. But even if Pgf = f, that is / is outside the loss subspace Hb, there 
may still be losses since all system degrees of freedom can be coupled. The analysis of the 
stored and dissipated energies, in view of the relations P6|) - (!5T]) . depends on the admittance 
operator 21 (cu). To study the properties of the admittance operator 21 (cu ) defined by ( 148]) 
we consider the block form (ITT)) and ([IB]) and represent cul — A (0) as 2 x 2 block operator 
matrix 



cul — A (P) = 
E 2 (ojJ) :=cuI 2 -(n 2 -il3B 2 ) 



S 2 (cu,P) 

-e* 



-e 

Si (cu) 
Sx(o;) 



(52) 



ojI\ — Qi 



where I 2 and Ji denote the identity operators on the spaces H B and H^, respectively. With 
respect to this block representation, the Schur complement of S 2 (cu, (3) in cul — A ((3) is 
defined as the operator 



S 2 (cu,P) = E 1 (cu)-e*E 2 (cu^y'Q, 



(53) 



whenever S 2 (cu, (3) is invertible. 

In what follows we assume the frequency cu ^ is not one of the resonance frequencies, 
that is cu 7^ pj, Nb + 1 < j < N. Then we know by Proposition [21] that the operators 
Si (a;), H 2 (cu,(3), S 2 (cu,f3), and cul — A ((3) are invertible for (3 3> 1. To simplify lengthy 
expressions we will suppress the symbols cu, (3 appearing as arguments in operators Hi (cu), 
E 2 (cu,(3), S 2 (cu,(3). Furthermore, the explicit formula based on the Schur complement is 
derived for the admittance operator 



21 (cu) = i[ul - A(f3)}~ 1 = i 



"2 + "2 ^^2 ®*^<2 



s 2 l 



(54) 
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A perturbation analysis at = oo of the admittance operator 21 (u), the results of which are 
summarized in Proposition l2~ij yields the following asymptotic expansion for /3 — > oo 



21 M 





o i-r 1 



+ w^p- 1 + o (/r 2 ) , w^^o, 



(55) 



where 



B 2 l > 0. 



(56) 



These asymptotics for the admittance operator lead to asymptotic formulas as (3 — > oo for 
the energy U, the dissipation power W^, and the quality Q factor which depend on whether 
Pgf = or Pgf 7^ 0. Namely, if P B f = 0, that is if / is inside the loss subspace H B , then 
Theorem (231 tells us that 



u = \{f, W + b^q (~^y~^e*B 2 i] /) /r 2 + o (/r 3 ) 

w dis = (f,B 2 i f) ^ + o(r 2 ), 

§(/, [5 2 - 2 + 5 2 - 1 e(s r 1 )*s r 1 0^ 2 - 1 ] /) 



(57) 



Q 



\CJ\ 



(/,^7) 



/r 1 + o (/r 2 ) 



and the leading order terms of U, Wdi S and Q are positive numbers. In particular, the quality 
factor Q — > as f3 — > oo. 

If p£f ^ then Theorem El tells us that 



u= 1 -(^p^f^p^f) + o(r l ) 



(58) 



(lis 



and the leading order term of U and is a positive and a nonnegative number, respectively. 
Furthermore, the quality factor is either infinite for (3 3> 1 (the case Wdi S = 0) or Q — >■ oo as 
/3 — t- oo. Moreover, 







-/3 + 0(l) 



(59) 



provided (/, W^~^ f) ^ 0, in which case the leading order term for Q is a positive number. 

Therefore the quality factor Q satisfies lim^oo Q = oo provided / has a non-zero projec 
tion on the no- loss subspace = H Q H B , whereas otherwise lim^oo Q = 0. 



3 An electric circuit example 

One of the important applications of our methods described above is electric circuits and 
networks involving resistors representing losses. A general study of electric networks with 
losses can be carried out with the help of the Lagrangian approach, and that systematic study 
is left for another publication. For Lagrangian treatment of electric networks and circuits we 
refer to [Gantl Sec. 9], [Goldl Sec. 2.5], [Pari] . 

We illustrate the idea and give a flavor of the efficiency of our methods by considering 
below a rather simple example of an electric circuit as in Fig. [TJ This example will show the 
essential features of two component systems incorporating high-loss and lossless components. 
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Li 



A/W 



R 



I, 




Figure 1: An electric circuit involving three capacitances Ci, C2, C12, two inductances Li, 
L 2 , a resistor i? 2 , and two sources E x , E 2 . 



To derive evolution equations for the electric circuit in Fig. [T] we use a general method 
for constructing Lagrangians for circuits, |Gantj Sec. 9], that yields 



rj, E\ . 2 L 2 . 2 

T = — 9i + y<? 2 > 



C7 



2Ci 



2C 



12 



2C< 



D R<1 -2 



(60) 



where T and {7 are respectively the kinetic and the potential energies, T—U is the Lagrangian, 
and R is the Rayleigh dissipative function. Notice that I\ = q± and I 2 = <j2 are the currents. 
The general Euler-Lagrange equations are, |Gantl Sec. 8], 



d dT dT 
dt dq\j dqj 



dU_ OR 

d qj d qj ' 



and for (1601) the general Euler-Lagrange equations therefore have the following form 



d_ 
dt 



Lxqi 



d_ 

dt 



L 2 q 2 



a 



■92 



Qi 



c 



c 



7i - 42) + fi, 



12 



12 



If we introduce 



Q 



L x 
L 2 

the system (IB"2"j) can be recast into 



Qi 

Q'2 

Ci 1 H 



Zi - 92) - ^2^2 + /a- 



/1 

/a 



°12 



°12 



°12 



°12 



R 




R 2 



Ld 2 t Q + R<9*Q + GQ = F. 
To provide for an efficient spectral study of the vector equation (IB"5l) we notice that 

L > 0, G > 0, 



(61) 



(62) 



(63) 
(64) 

(65) 
(66) 
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L^RLri 



I_-*GL - * 



and introduce 

Q = iJQ, F = [r*F, R 
The vector equation (!65|) is transformed then into 

d 2 t Q + Rd t Q + GQ = F. 
To rewrite the above second-order ODE as the first-order ODE we set 



X 



allowing one to recast the vector equation (l68l as 



y 




' d t Q ' 


X 







that is y = d t Q, 



x 



d t X 



-R 



F 




where 













-$11 


-$12 


— R -$ " 







—R2L2 1 


-$12 


-$22 


$ 




$11 


$12 












$12 


$22 









Consequently, the general Euler-Lagrange equations fl62|) for the electric circuit in 
transformed into the canonical form (jSj), namely 



(67) 
(68) 

(69) 
(70) 

(71) 
Fig. [T]are 



d t X = -i (O - i/35) X + 



F 




with the system operator 



n 



A{fi) = Q-if3B, 0>O, 






_i$ " 


, $ = 


" $ n $ 12 " 


i$ 





_ $12 $22 _ 



F 










r- 1 







and (3 



R 2 t 



> 0, 



where r > is a unit of time. 



(72) 
(73) 

(74) 



As we can see this electric circuit example fits within the framework of our model. Indeed, 
since resistors represent losses, this two component system consists of a lossy component and 
a lossless component - the right and left circuits in Fig. (TJ respectively. 

The next two sections are devoted to the analysis of this electric circuit both theoretically 
and numerically in the high-loss regime 3> 1 using the methods developed in this paper. 
For this purpose the following properties of the matrix $ in ( 1731) are useful: 



$ = ^Tr($ 2 ) + 2v/det ($ 2 )^) (y/det ($ 2 )/ 2 + $' 



(75) 



Li 



C\ C12 / 



1 



$11, $22 > 0, 
where a/ 7 denotes the positive square root. 



/L1VL2C12 



L2 \ C2 

$12 e 



C12 ) 
^\{0}, 



> 0, 
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3.1 Spectral analysis in the high- loss regime 

In this section, a spectral analysis of the electric circuit example in Fig. [I] in the high-loss 
regime is given using the main results of this paper. 



3.1.1 Perturbation analysis 

The finite dimensional Hilbert space is H = C 4 under the standard inner product (•, • 
is decomposed into the direct sum of invariant subspace of the operator B > in (1T4T) , 



It 



H = H 



B 



H 



where Hb = rani?, H B = ker B are the loss subspace and no-loss subspace with dimensions 
Nb = 1, N — Nb = 3 and orthogonal projections 



P, 



B 





10 







p 



B 



10 



10 

1 



respectively. The operators Q and B, with respect to this direct sum, are 2x2 block operator 
matrices 

n 



' n 2 


e " 


, B = 


' B 2 


" 


e* 











where Q 2 := Pb^Pb\h b '■ Hb — > Hb and B 2 := PbBPb\ Hb : Hb —> Hb are restrictions of the 
operators f2 and B, respectively, to loss subspace Hb whereas f2i := PgClPM ± : H^ — y Hg 

is the restriction of Q to complementary subspace Hg. Also, G : H^ —y Hb is the operator 
:= P B Q,Pk\ ± whose adjoint is given by 6* = PkQ,P B \ : H B -> H B . Moreover, 

according to our perturbation theory the operator Q*B 2 ~ 1 Q : H B — >■ H B plays a key role in 
the analysis. These operators act on the 4x1 column vectors in their respective domains as 
matrix multiplication by the 4x4 matrices 



n 2 = 0, B 2 = B, fix 



e 



B? 













r 









i$i 2 










-r$ n 


_i$ 














i$n 























12 





"i$22 






e* 



o 










i$i2 

i$22 












T$i2$22 














T$12$22 



where by " =" we mean equality as functions from the domain of the operator on the LHS of 
the equal sign. 
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The operators fii and B 2 in this example have only simple eigenvalues and so we will 
use Corollary [12j We introduce below a fixed orthonormal basis {wj}^ =1 diagonalizing the 
operators Q\ and B 2 and then determine the values Q, pj, dj from the relations 

B 2 Wj = CjWj'j Pj = ( w j, &>Wj) for j = 1; 

Qiibj = Pjihj, dj = (wj, 0* B 2 Qwj] for 2 < j < 4. 



In particular, 




1 





Ci = r 



-1 



Pi = 0; 



(76) 



^2 






-$ 12 



r($f 2 -$ n $ 22 ) 

^ = ' * = $ 2 + $ 2 > > 

1 1 1 T ^12 



(77) 



w 3 



V2 y/$ 2 u + $ 



2 

12 



-iV^ii + * 


$11 
$12 



2 

12 



P3 



r$f 2 ($11 + $22) 



W4 = w 3 , p 4 = — p 3 < 0, d/i = d 3 > 0. 



By Theorem [5] and Corollary [T2l of this paper it follows that in the high-loss regime ^> 1, 
the system operator A ((3) = Q — i(3B is diagonalizable and there exists a complete set of 
eigenvalues and eigenvectors satisfying 

A(f3)w J (f3) = Q(f3)w j (f3), 1 < j < 4, 

which splits into two classes 

high-loss: Q (/3) , ^ (/3) , j — 1; (78) 
low-loss: Ci , ™j (0) , 2 < j < 4, 

having the following properties. 

High-loss modes. The high-loss eigenvalue has a pole at /3 = 00 whereas its eigenvector 
is analytic at = 00, having the asymptotic expansion 

Ci(/3) = -iCi/3 + Pi + 0(r 1 ), 6 > 0, Pl eR, (79) 

w) = ^1 + (r 1 ) , 

as (3 00. The vector w\ is an orthonormal basis of the loss subspace Hb and 

Bwi = (iwi, pi = (wi, ttwi) . 

Low-loss modes. The low-loss eigenvalues and eigenvectors are analytic at = 00, 
having the asymptotic expansions 



08) = p 3 - i^/T 1 + O (/T 2 ) , p 3 e R, d 3 - > 0, 
™i (/?) = Wj + u^/T 1 + O (/T 2 ) , 2 < j < 4, 
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as (3 — > oo. The vectors Wj, 2 < j < 4 form an orthonormal basis of the no- loss subspace Hg 
and 

Bvjj = 0, pj = (wj, ttuij) , dj = (w^, Bw { r l)Ss j , 2 < j < 4. 
3.1.2 Overdamping and symmetries of the spectrum 

The phenomenon of overdamping (also called heavy damping) is best known for a simple 
damped oscillator. Namely, when the damping exceeds certain critical value all oscillations 
cease entirely, see, for instance, |Pain[ Sec. 2]. In other words, if the damped oscillations 
are described by the exponential function e -1 *'* with a complex constant ( then in the case 
of overdamping (heavy damping) Re ( = 0. Our interest in overdamping is motivated by the 
fact that if an eigenmode becomes overdamped then it will not resonate at any frequency. 
Consequently, the contribution of such a mode to losses becomes minimal, and that pro- 
vides a mechanism for the absorption suppression for systems composed of lossy and lossless 
components. 

The treatment of overdamping for systems with many degrees of freedom involves a num- 
ber of subtleties particularly in our case when the both lossy and lossless degrees of freedom 
are present. We have reasons to believe though that any Lagrangian system with losses ac- 
counted by the Rayleigh dissipation function can have all high-loss eigenmodes overdamped 
for a sufficiently large value of the loss parameter (3. In order to give valuable insights into 
far more general systems, we focus on the electric circuit example in Fig. 1 giving statements 
and providing arguments on the spectral symmetry and overdamping for the circuit. This 
analysis is used in the next section to interpret the behavior of the eigenvalues of the circuit 
operator A ((3). 

Our first principal statement is on a symmetry of the spectrum of the system operator 
A (/3) with respect to the imaginary axis. 

Proposition 1 (spectral symmetry) Let A {(3) denote the system operator for the 
electric circuit given in Fig. [H Then for each (3 > 0, its spectrum o (A (/?)) lies in the lower 
half of the complex plane and is symmetric with respect to the imaginary axis, that is 

a (A (f3)) = -a (A (/?)). (81) 

Moreover, except for a finite set of values of (3, the system operator A ((3) is diagonalizable 
with four nondegenerate eigenvalues. 

Proof. From the asymptotic analysis in (J79l) -( J771) it follows that all the eigenvalues of 
A (f3) must be distinct for ^> 1. Now the operator Q — i(3B, (3 e C is analytic on C and so, 
by a well-known fact from perturbation theory |Bau85t p. 25, Theorem 3; p. 225, Theorem 
1] , its Jordan structure is invariant except on a set S C C which is closed and isolated. These 
facts imply the system operator A ((3) is diagonalizable with distinct eigenvalues except on 
the closed and isolated set S fl [0, oo) which must be bounded since the eigenvalues of A ((3) 
are distinct for [3 ^> 1. In particular, this implies S PI [0, oo) is a finite set. This proves that 
except for a finite set of values of f3, the system operator A {(3) is diagonalizable with four 
nondegenerate eigenvalues. 

Next, since A ((3) = Q — i(3B in (172|) is a system operator satisfying the power dissipation 
condition B > then it follows from Lemma [271 in Appendix E] that if (3 > then Im^ < if 
( is an eigenvalue of A ((3), i.e., the spectrum a (A ((3)) lies in the lower half of the complex 
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plane. Finally, one can show that det ((I - A (/?)) = det L" 1 det (C 2 L + (iR - G) for all (eC, 
(3 > 0. Moreover, from our assumptions /3 > 0, r > and L, G > it follows that the 2x2 
matrices L, R, and G must have real entries. Using these two facts we conclude 

det ((I - A (0)) = det L" 1 det (( 2 L + (iR - G) = (82) 

= det L" 1 det ((-C) 2 L + (-()iR - g) = det ((-£)/- A (/?)), 

and, hence, (IHTj) holds. ■ 

Corollary 2 (eigenvalue symmetry) Lei X fre any open interval in (0, oo) with the prop- 
erty that the eigenvalues of system operator A ((3) are nondegenerate for every (3 G X. Then 
there exists a unique set of functions Q : X — > C, j = 1,2,3,4 which are analytic at each 
(3 G X and whose values Ci (P), C2 ifi), C3 (P)> C4 (£0 ore the eigenvalues of the system operator 
A(j3). Moreover, there exists a unique permutation k : {1,2,3,4} 1 — > {1,2,3,4} depending 
only on the interval X such that for each j = 1, 2, 3, 4, 

(3 (P) = -C«y) (#> for every G X. (83) 

Proof. It is a well-known fact from perturbation theory for matrices depending analyti- 
cally on a parameter |Bau85j . that simple eigenvalues can be chosen to be analytic locally in 
the perturbation parameter and analytically continued as eigenvalues along any path in the 
domain of analyticity of the matrix function which does not intersect a closed and isolated 
set of singularities. These singularities are necessarily contained in the set of parameters in 
the domain where the value of matrix function has repeated eigenvalues. The proof of the 
first part of this corollary now follows immediately from this fact and the fact the high-loss 
and low-loss eigenvalues of A (/3) are meromorphic and analytic at (3 = 00 , respectively, with 
distinct values for f3 3> 1. The existence and uniqueness of the permutation is now obvi- 
ous from this and symmetry of the spectrum described in the previous proposition. This 
completes the proof. ■ 

Corollary 3 (overdamping) Let Q ((3), j = 1,2,3,4 be the high-loss and low-loss eigen- 
values of the system operator A ((3) given by J7^)-(f7^). Then, in the high-loss regime (3^1, 
these eigenvalues lie in the lower open half-plane and, moreover, the eigenvalues (j (j3), 
j — 1,2 are on the imaginary axis whereas the eigenvalues (j ((3), j = 3,4 lie off this axis 
and symmetric to it, i.e., (4 (0) = —( 3 (j3). 

Proof. First, it follows the asymptotic analysis in (!79|) -( !77|) that there exists a f3 > 
such that Q ((3), j = 1,2,3,4 are all the eigenvalues of system operator A (f3) and are 
distinct for every (3 G (f3 , 00). By the previous corollary there exists a unique permutation 
k : {1,2,3,4} 1 — > {1,2,3,4} depending only on the interval (/3 , 00) such that for each 
j = 1, 2, 3, 4, the identity flM} for every (3 G (/3 , 00). 

Next, we will now show that for this permutation we have k(1) = 1, k(2) = 2, k(3) = 4, 
«; (4) = 3. Well, consider the the asymptotic expansions of the imaginary and real parts 
of high-loss and low-loss eigenvalues. First, lim^oo Im £1 (j3) = —00, lim^oo Im Q (/3) = 
0, j = 2,3,4 and since CiW) = — Ck(i) {P) these properties imply k(1) = 1. Second, 
lim^oo Re C2 (P) = 0, hm ^^ Re C 4 (P) = p4 = ~p3 = ~ hm^^ Re C 3 (P) with p 3 > and 
since Cj (P) = ~ Gc(j) (P) with n(j) 7^ 1 for j = 2,3,4 these properties imply re(2) = 2, 
k(3) = 4, k(4) = 3. 
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To complete the proof we notice that since k(1) = 1, k(2) = 2, k(3) = 4, k (4) = 3 then 
for (3 3> 1 we have Q (f3) = —Q (f3), j = 1,2 and C4 (P) = —(3 (/?)• The proof now follows 
immediately from this and the facts — Re£ 4 (/3) = Re( 3 (f3) = p 3 + 0(/3™ 2 ), Im£ 4 (/3) = 
ImC 3 = -d 3 f3~ l + O (/3 -3 ) as f3 -)■ 00 where p 3 , d 3 > 0. ■ 

Remark 4 TTie boundary of the overdamping regime known as critical damping, corresponds 
to a value of the loss parameter (3 = (3 > at which the system operator A ((3) develops a 
purely imaginary but degenerate eigenvalue Co- The spectral perturbation analysis of A ((3) 
in a neighborhood of the point (3 = f3 is theoretically and computationally a difficult problem 
since it is a perturbation of the non- self- adjoint operator A (f3o) with a degenerate eigenvalue. 
This type of perturbation problem was considered in J Welti/ where asymptotic expansions of the 
perturbed eigenvalues and eigenvectors were given and explicit recursive formulas to compute 
their series expansions were found lWelt\ Theorem 3.1], under a generic condition I Welt\ p. 
2, (1.1)]. In particular, this condition is satisfied for the system operator A (f3) at the point 
— A) f or the degenerate eigenvalue Co since 

^ det (CI - A (/?)) | (Cj/3 ) =(CoA ) = ir- 1 ^ 3 - ir- 1 + Co ^ 0. 
3.2 Numerical analysis 

In order to illustrate the behavior of the eigenvalues of the system operator for the circuit in 
Fig. [T]we fix positive values for capacitance Ci, C 2 , Cu, inductances Li, L 2 , and the unit 
of time r. Once these are fixed, the system operator A{j3) is computed using fl64|) . (|67|) . and 
(|72l) -(j74 l . These values constrain the magnitude of the resistance R 2 of the corresponding 
circuit in Fig. [1] to be proportional to the dimensionless loss parameter (3 since it follows 
from (El) that 

R 2 = —P. (84) 

r 

The high-loss regime f3 3> 1 is associated with the right circuit in Fig. [1] experiencing huge 
losses due to the resistance R 2 ^> 1 while the left circuit remains lossless. In particular, 
each choice of these values provides a numerical example of a physical model with a two 
component system composed of a high-loss and lossless components. 
For the numerical analysis in this section we chose 

Ci := 2, C 2 := 3, C 12 := 4, L x := 5, L 2 := 6, r := 1. (85) 

All graphs were plotted in Maple® using these fixed values and with the loss parameter in 
the domain < (3 < 10. 



On Figure [21 In Fig. |2] is a series of plots which compares the imaginary part of each of 
the eigenvalues Q ((3), 1 < j < 4 of the system operator A((3) for the electric circuit in Fig. 
[T]to the quality factor Q of their corresponding eigenmodes. To plot the quality factor as 
a function of the loss parameter (3 we have used formula (I3!)]) . As is evident by this figure, 
there is clearly modal dichotomy caused by dissipation. 

Indeed, for the high- loss eigenpair Ci ((3), w\ (/?) we can see in plots (a)-(b) that as the 
loss parameter (3 grows large so too does Im Ci (/?) whereas the quality factor Q[wi (/?)] = 
— I t^t§y of the mode goes to zero as predicted by our general theory (cf. (|99|) of Proposition 
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Figure 2: (a)-(d) For each high- loss eigenvalue Q (j3), j = 1 and low- loss eigenvalue Q (f3), 
2 < j < 4 of the system operator A((3) for the electric circuit in Fig. [I] with the values 
(1B4|) and fl85l) . comparing its imaginary part Im(j ((3) to the quality factor Qj = — | 
of any one of its eigenmodes (not shown is Qj — )■ oo as /3 — > 0). For /3 ^ (3 w 0.57282 
(critical damping) and j = 1,2, the eigenmodes with eigenvalue (/?) are overdamped since 
Re£j ((3) = 0. The overdamping phenomenon is manifested graphically in (a), (c) with the 
cusps in the curves at the intersection of vertical dotted line (3 = (3$ and (b), (d) with the 
curves on the line Q = for (3 ^ (3$. (a) Im£i (/3) vs. (3. (b) Q\ vs. (3. Due to overdamping, 
Qi = for (3 > 0o. (c) IraQ (/?) vs. (3, j = 2,3,4. The curves Im£ 2 (/3) and Im£ 4 ((3) are 
shown in blue and green, respectively. Due to eigenvalue symmetry, the curves Im^ ((3) and 
Im£ 4 ((3) cannot be distinguished in this plot and Im£ 2 ((3) = Im£i (/3) for (3 < (3q. (d) Qj 
vs. (3, j = 2, 3, 4. The curves Q = Q2 and Q = Q4 are shown in blue and green, respectively, 
and because of overdamping Q 2 = for (3 ^ (3 Q . Due to eigenvalue symmetry, the curves 
Q = Q3 and Q = Q4, cannot be distinguished in this plot. 
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Figure 3: (a)-(d) A comparison of the real and imaginary parts of the low-loss eigenvalue 
£ 3 (/3) to its truncated asymptotic expansion ( 3 = p 3 — i<i 3 /3 _1 for the values p 3 , d 3 predicted 
by our theory As evident from these plots, £3 (f3) ~ ( 3 ((3) for (3 large, (a) Re£ 3 (/?) vs. /3. 
(b) Re£ 3 (/3) vs. (3 (the solid line) and Re (3 vs. /3 (the diagonal crosses), (c) Im£ 3 (/3) 
vs. /3. (d) Im£ 3 (j3) vs. /3 (the solid line) and Im£ 3 (/3) vs. j3 (the diagonal crosses). 
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Hand (11 101) of Proposition [T4"j) . In fact, by our results on the overdamping phenomenon 
described in Corollary E] we know that it is exactly zero for all (3 > (3 , where (3 denotes 
the boundary of the overdamped regime as discussed in Remark [U For the fixed values in 
(185)) . O ~ 0.57282 and we have place a vertical dotted line in each of the plots in Fig. [2] to 
indicate this boundary. 

The behavior is quite different for the low-loss eigenpairs Q ((3), Wj ((3), 2 < j < 4. The 
plots (c)-(d) show that as the loss parameter (3 grows large, the values Im Q j = 2,3,4 
all become small with the quality factors Q[w 3 (/?)] = — | ^^[^ and Q[u; 4 ((3)] = —§ ^^03) 
of the eigenmodes w 3 (f3) and u> 4 ((3) becoming large as predicted by our theory (cf. (1991) of 
Proposition [71 ( I112p of Proposition HH and formulas (177j) ). The quality factor Q[w 2 (/?)] = 
~\ ^!mC2(^ °^ ^ e l° w "l° ss m °de w 2 ((3) becomes zero for (3 > (3 , again a fact which is 
predicted for this electric circuit from the overdamping phenomenon described in Corollary 

El 



On Figure [31 Figure El compares the low-loss eigenvalue £3 {(3) of the system operator A {[3) 
to the truncation ( 3 of its asymptotic expansion as predicted by our theory in (THUj) and 
( 1771) . namely, 

C3 (/?) « C3 W) = P3~ = J&n + " P » 1- 

^11 1 ^12 

The plots (a) and (c) in the figure are the real and imaginary parts, respectively, of the 
eigenvalue (3 and plots (b) and (d) are the real and imaginary parts, respectively, of both 
the eigenvalue (3 ((3) and the truncation of its asymptotic expansion (3 ((3). 



On Figure [4] In Fig. H]we have the images in the complex plane of the eigenvalues (j 
1 < j < 4. This figure displays the spectral symmetry and overdamping phenomena as 
predicted in Proposition [1] and Corollaries El El According to Lemma |27J and as evident 
in the figure, the eigenvalues lie in the lower half-plane. By Theorem [151 these eigenvalues 
converge to the real axis as (3 — > and, in particular, to the eigenvalues of the frequency 
operator Q. 

In plot (a) we see the images of the eigenvalues Q ((3), j = 1,2 in the complex plane. 
We observe that overdamping does occur but only for these two eigenvalues. Indeed, as the 
loss parameter (3 increase from zero these two eigenvalues eventually merge on the negative 
imaginary axis when (3 = /?o ~ 0.57282 and stay on this axis for all (3 > (3o with Im£i {(3) — > 
— 00 and Im(2 (/?) — > as (3 — )• 00. 

In plot (b) we see the images of the eigenvalues Q ((3), j = 3, 4 in the complex plane. This 
plot shows the eigenvalue symmetry for the system operator A ((3) for the electric circuit which 
forces these eigenvalues to satisfy £4 = — (3 to lie off the imaginary axis, and to be in 
the lower open half-plane for all (3 > with hn.Q ((3) — > as (3 — )• 00 for j = 3, 4. 
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Figure 4: (a)-(b) Image in the complex plane of the eigenvalues Q (ft), 1 < j < 4 of the 
system operator A(/3) for the electric circuit in Fig. [I] with the values ( 1841 and ( )85l) . The 
power dissipation condition implies lm(j(ft) < for ft > 0, as evident in the figure, (a) 
Image of the high-loss eigenvalue Ci (ft) and low-loss eigenvalue £2 (ft) displayed in red and 
blue, respectively. For the purpose of comparison, the view is restricted to a box around 
the image of (2 (ft)- Off the imaginary axis the blue curve is symmetric about this axis to 
the red curve due to the eigenvalue symmetry ( 2 (ft) = — Ci (ft) f° r ft < fto ~ 0.57282. The 
two curves intersect, for ft = fto, on the negative imaginary axis and due to overdamping 
stay there for all ft ^ fto- Moreover, Im(i (ft) — > —00 and Im£ 2 (ft) — > as ft — > 00. (b) 
Image of the low-loss eigenvalues ( 3 (ft) and (4 (ft) displayed in orange and green, respectively. 
The green curve is symmetric about the imaginary axis to the orange curve due to the 
eigenvalue symmetry £ 4 (ft) = —( 3 (ft) for all ft ^ 0. Moreover, Im£ 4 (ft) = Im£ 3 (ft) — > 
and - Re ( 4 (ft) = Re C3 (ft) -> p3 as ft ->• 00, where p 3 w 0.40825. 
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4 Perturbation analysis of the system operator 



This and the following sections are devoted to a rigourous perturbation analysis of the eigen- 
values and eigenvectors of the system operator from (jSj), 

A (/3) — i/3B, 0>O, (86) 

in both the high-loss regime, /3 > 1, and the low-loss regime, -C 1. We mainly focus on the 
high-loss regime and our goal is to develop a mathematical framework based on perturbation 
theory for an asymptotic analytic description of the effects dissipation have on the system dHJ) 
including modal dichotomy, i.e., splitting of eigenmodes into two distinct classes according to 
their dissipative properties: high-loss and low-loss modes. This framework and its rigorous 
analysis provides for insights into the mechanism of losses in composite systems and in ways 
to achieve significant absorption suppression. 

The rest of this section is organized as follows. We first recall the basic assumptions, 
definitions, and notations from earlier in this paper regarding the system operator ( 186]) . the 
quality factor Q and power of energy dissipation Wdi B associated with its modes. In the 
next section we state our main results on the perturbation analysis of the eigenvalues and 
eigenvectors for this operator A(/3). The result for the high-loss regime, ^> 1, and the 
low-loss regime, /? <C 1, are placed in separate sections. Finally, we prove the statement of 
our main results in Section [6] 

The system operator A(/3) in (|86|) . for each value of the loss parameter (3, is a linear 
operator on the finite dimensional Hilbert space H with N := dimP and scalar product 
(•,•). The frequency operator Q and the operator associated with dissipation B are self- 
adjoint operators on H. The operator B satisfies the power dissipation condition (Q and 
the loss fraction condition (|T5|) . namely, 

B > 0, < 5 B < 1 (87) 

where Nb ■— rankP denotes the rank of the operator B and 5b '■= -# is referred to as the 
fraction of high-loss modes. 

The range of the operator B, i.e., the loss subspace, is denoted by Hb and the orthogonal 
projection onto this space is denoted by Pg. It follows immediately from these definitions 
and the fact B is self-adjoint that 

H = Hb © Hg (88) 

where H^, i.e., the no-loss subspace, is the orthogonal complement of Hb in H and is the 
kernel of B with the orthogonal projection onto this space given by Pg :— I — Pg. In 
particular, 

H B = ranP, N B = dimH B , (89) 
H% = ker P, N - N B = dim H%. 

The energy U [w] , power of energy dissipation W<n s [w] , and the quality factor Q [w] of an 
eigenvector w of the system operator A (/3) with eigenvalue ( is 

I - (yj yj) 

U[w] = -(w,w), W dis [w} = (w,(3Bw), Q[w] = \Re(\ aV ' \ , (90) 
z [w, prSw) 
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where Q [w] is said to be finite if Way [w] ^ 0. In Appendix [8] we show that 

ImC = - ( 7^w W dis [w] = -2lmCU[w], QH^-JtT' ( 91 ) 
(w,w) 2 Im<; 

where Q [«;] is finite if and only if Im £ ^ 0. 
4.1 The high-loss regime 

We begin this section with our results on the perturbation analysis of the eigenvalues and 
eigenvectors for this operator A ((3) for the high- loss regime in which ^> 1. 

Theorem 5 (eigenmodes dichotomy) Let Q, 1 < j < Nb be an indexing of all the 
nonzero eigenvalues of B (counting multiplicities) where Nb = rank B. Then for the high- 
loss regime (3^1, the system operator A((3) is diagonalizable and there exists a complete set 
of eigenvalues Q ((3) and eigenvectors Wj ((3) satisfying 

A ((3) Wj ((3) = Q ((3) Wj ((3), l<j< N, (92) 

which split into two distinct classes of eigenpairs 

high-loss: Q (P) , Wj (/9) , l<j< N B ; (93) 
low-loss: Q {(3) , Wj (13), N B + 1 <j < N, 

having the following properties: 

(i) The high-loss eigenvalues have poles at (3 = oo whereas their eigenvectors are analytic 
at (3 = oo. These eigenpairs have the asymptotic expansions 

( 3 ((3) = -i( J f3 + p j + 0((3- 1 ) } 4>0, (94) 
Wj ((3) = Wj + (/T 1 ) , 1 < j < N B 

as (3 oo. The vectors Wj, 1 < j < N B form an orthonormal basis of the loss subspace 
Hb and 

Biij = (jWj, pj = (vjj,Qwj) , (95) 

for 1 < j < Nb- 

(ii) The low-loss eigenpairs are analytic at (3 = oo and have the asymptotic expansions 

Q ((3) = P j - idj/3- 1 + O (/T 2 ) , Pj e R, dj > 0, (96) 
Wj {(3) = Wj + w^P' 1 + O (/T 2 ) , N B + l<j<N 

as (3 — > oo. The vectors Wj, Nb +1 < j < N form an orthonormal basis of the no-loss 
subspace and 

Bwj = 0, pj = (wj, ttwj) , dj = Bwf l) ^j (97) 

for N B + 1 < j < N. 
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Corollary 6 (eigenmode expulsion) The projections of the eigenvectors Wj ((3), 1 < j < 
N onto the loss subspace H B and the no-loss subspace have the asymptotic expansions 

P BW] (?) = Wj + O (/T 1 ) , P&Wj (?) = O (/T 1 ) , l<j<iV fl ; (98) 
Pjtwj ifi) = Wj + (/T 1 ) , P B w 3 {fi) = O (I3- 1 ) , N B + l<j<N 

as (3 — ^ oo. 

Proposition 7 (eigenfrequency expansions) The functions KeQ ((3), l<j<Nat/3 = 
oo are analytic and their series expansions contain only even powers of (3~ l . The functions 
Im(j ((3), l<j<Nat(3 = oc have poles for 1 < j < N B , are analytic for N B + 1 < j < 
N, and their series expansions contain only odd powers of f3~ l . Moreover, they have the 
asymptotic expansions 

Re( 3 ((3)=p 3 +0(f3- 2 ), lm( J (f3) = -Q(3 + 0((3- 1 ), 1 < j < N B ; (99) 
ReQ(f3)=p 3 +0(f3- 2 ), ImQ (f3) = -d^' 1 + O (f3~ 3 ) , N B + l<j<N 
as (3 —7- oo. 

Proposition 8 For each j = 1,...,N and in the high-loss regime > 1, the following 
statements are true: 

1- Ifl<j< N B then ImQ (f3) < 0. 

2. If N B + 1 < j < N then either ImQ ((3) = or ImQ (J))<0. Moreover, ImQ (0) = 
if and ifQ (f3) = pj. 

3. If N B + 1 < j < N then Wj $ ker (pjl - Q) if and only if d 3 ^ 0. 

4- If N B + 1 < j < N then Wj ^ ker Q if and only if pj ^ or d 3 • ^ 0. 

Remark 9 Typically, one can expect that the asymptotic expansion of the low-loss eigenval- 
ues have dj ^ 0, for j = N B + 1, . . . , N. Indeed, if this were not the case then Theorem 
and the previous proposition tell us that the intersection of one of the eigenspaces of the 
operator Q with the kernel of the operator B would contain a nonzero vector. And this is 
obviously atypical behavior. 

Remark 10 An important subspace which arises in studies of open systems in \Liv\ pp. 
27-28] as well as in \FigShi , Sec. 4-1] is 

O n (ran B) = Sp&n{Q n Bu : u e H, n = 0, 1, . . .} , (100) 

where it is called the orbit and is the smallest subspace of H containing rani? that is invariant 
under Q. It is shown there that the orthogonal complement Oq (ranB) 1 ' = HQ On (rani?) = 
fl n >o ker (BQ n ) is the subspace which is invariant with respect to the operator Q — i(3B and 
the restriction Q — i(3B\ H ± is self-adjoint. Consequently, the evolution over this subspace in 
entirely decoupled from the operator B and there is no energy dissipation there. Moreover, 
the importance of the orbit to the perturbation analysis is that in applications it's often not 
hard to see that 

O n (™nB) = H. (101) 

in this case, by Remark^ we know that the asymptotic expansion of the low-loss eigenvalues 
have dj ^ 0, for j = N B + 1, . . . , N. 
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For computational purposes the next proposition and its corollary are important results. 
We first recall some notation. The orthogonal projection onto the loss subspace ker B = Hb 
is Pb and is the orthogonal projection onto the no-loss subspace rani? = Hg. Then from 
(US) we will need the operators B 2 = P b BP b \ Hb : H B -> H B ,Sl\= P B ^P B \ H± '■ h b ~* h bi 
and 6 = P B QPM„ ± : H% -> H B whose adjoint is 9* = P£QP B \„ :H b ^H^. 

Proposition 11 (asymptotic spectrum) The following statements are true: 

1. In the asymptotic expansions (OM) for the high-loss eigenpairs, the coefficients Q, Wj, 
1 < J < Nb form a complete set of eigenvalues and orthonormal eigenvectors for the 
operator B 2 with 

B 2Wj = Qwj, l<j< N B . (102) 
In particular, B 2 is a positive definite operator as is its inverse B 2 X , i.e., 

B 2 > 0, B 2 X > 0. (103) 

2. In the asymptotic expansions $UbT) for the low-loss eigenpairs, the coefficients pj, Wj, 
N B + 1 < j < N form a complete set of eigenvalues and orthonormal eigenvectors for 
the self-adjoint operator Q\ with 

SliWj = pjWj, N B + 1 < j < N. (104) 

3. The coefficients dj, Nb + 1 < j < N in the asymptotic expansions of the low-loss 
eigenvalues are given by the formulas 

dj = (wj, e*B 2 l ewj) , N B + 1 < j < N. (105) 

Corollary 12 (computing expansions) The following statements give sufficient condi- 
tions that allow computation ofvjj, Q, pj, and dj in the asymptotic expansions of the eigen- 
pairs: 

1. If the eigenvalues of B 2 are distinct and qj, 1 < j < Nb is any indexing of these 
eigenvalues then in Theorem^ after a possible reordering of the high-loss eigenpairs in 
< fP3|) . the coefficients in the asymptotic expansions (94\ ) are uniquely determined by the 
relations 

(j = 9, B 2 Wj = Qwj, \\wj\ \ = 1, pj = (wj, ttwj) , l<j< N B . 

2. If the eigenvalues of Q± are distinct and Qj, Nb + 1 < j < N is any indexing of these 
eigenvalues then in Theorem^ after a possible reordering of the low-loss eigenpairs in 
( ESP , the coefficients in the asymptotic expansions If9h}) are uniquely determined by the 
relations 

Pj = Qj, VliWj = PjWj, \\wj\\ = 1, dj = (vjj, 9* B^Qvjj) , N B + 1 < j < iV . 

The next two propositions give the asymptotic expansions as ft — > oo of the energy, power 
of energy dissipation, and quality factor for the high- loss and low- loss eigenvectors Wj (j3), 
l<j<N. 
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Proposition 13 (energy and dissipation) The energy for each of the high-loss and low- 
loss eigenvectors have the asymptotic expansions 

U[w j (P)] = ± + 0(/3- 1 ), l<j<N (106) 

as (3 — >■ oo. The power of energy dissipation for the high-loss and low-loss eigenvectors have 
the asymptotic expansions 

W dls [ Wj (/?)] = Q(3 + 0(1), 1 < j < N B ; (107) 
W dis [ Wj (/?)] = d,(3' 1 + O (/T 2 ) , N B + l<j<N 

as (3 oo. In particular, 

lim^M/))] = {~ t! llfit^ N . (108) 

Proposition 14 (quality factor) For each j — 1, . . . , N , the following statements are true 
regarding the quality factor of the high-loss and low-loss eigenvectors: 

1. For (3^1, the quality factor Q [wj (j3)] is finite if and only iflmQ ((3) ^ 0. 

2. If the quality factor Q [wj (/?)] is finite for (3^1 then it is either analytic at [3 = oo or 
has a pole, in either case its series expansion contains only odd powers of (3~ l and, in 
particular, 

lim Q [wj (/?)] = or oo. 

3. The quality factor of each high-loss eigenvector is finite for j3 >> 1 and has the asymp- 
totic expansion 

Q[wj(0)] = ^P- 1 + O{l3- s ), l<j<N B (109) 
as (3 — >■ oo. In particular, 

lim Q [ Wj (/?)] = 0, l<j<N B . (110) 

0— >oo 

4- If j € {N B + 1, . . . , N} and dj ^ then the quality factor of the low-loss eigenvector 
Wj ((3) is finite for (3 ^> I and has the asymptotic expansion 

QHm = \^+o{(3- 1 ) an) 



2 dj 



as (3 oo. In particular, 



lim QK = (112) 
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4.2 The low-loss regime 

We now give our results on the perturbation analysis of the eigenvalues and eigenvectors for 
the system operator A ((3) in the low-loss regime < j3 < 1. The focus of this paper is 
on the high-loss regime and so we do not try to give results as general as those in previous 
section. Instead, the goal of this section is to show the fundamentally different asymptotic 
behavior in the low-loss regime compared to that of the high-loss regime. 

Theorem 15 (low-loss asymptotics) Letujj, 1 < j < N be an indexing of all the eigenval- 
ues ofQ (counting multiplicities). Then for < /3 < 1, the system operator A ((3) = Q — i(3B 
is diagonalizable and there exists a complete set of eigenvalues (j ((3) and eigenvectors Vj (f3) 
of A (f3) satisfying 

A (P) Vj (/?) = Cj (P) v j (P) , l<j<N (113) 
with the following properties: 

(i) The eigenvalues and eigenvectors are analytic at (3 = and have the asymptotic expan- 
sions 

( j (/3) = uj j -io- j !3 + 0(l3 2 ), Wj-6 1, <Tj>0, (114) 
Vj (/3) =Uj + 0(/3), 1 < j < N 

as (3 — y 0. The vectors Uj, 1 < j < N form an orthonormal basis of eigenvectors of Q 
and 

fluj = ujjUj, Oj = (uj, Buj) , 1 < j < N. (H5) 

Corollary 16 (energy and dissipation) The energy and power of energy dissipation of 
these eigenvectors have the asymptotic expansions 

U [Vj (/?)] = 1 + (/T 1 ) , W dls [vj (/?)] = a 3 (3 + O (f3 2 ) , l<j<N (116) 
as f3 — > 0. In particular, 

lim W dls [vj (/?)] = 0, 1 < j < N. (117) 

Corollary 17 (quality factor) The quality factor of each of these eigenvectors has the 
asymptotic expansion 

Qhm = l—f3- l + 0(f3), (118) 
as (3 — > 0, provided o~j ^ 0, in which case it has the limiting behavior 

lim Q [vj (/3)] = I °? (119) 

/3->o J y o if Uj = o. v ' 

Remark 18 Typically, one can expect that the asymptotic expansion of these eigenvalues 
have Oj ^ 0, for 1 < j < N . Indeed, if this were not the case then it would follow from the 
assumption B > and U15\) of Theorem\T^ that the intersection of one of the eigenspaces of 
the operator Q with the kernel of the operator B would contain a nonzero vector. And this 
is obviously atypical behavior as mentioned previously in Remark^ 
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These results show that in the low- loss regime < f3 <C 1, all the modes behave as 
low-loss modes since the power of energy dissipation is small and typically the quality factor 
is very high. In contrast, the high- loss regime /3 > 1 has both a fraction < 5 b < 1 of 
high-loss modes and a fraction < 1 — 5b < 1 of low-loss modes. The behaviour of the 
low-loss modes in either regime is similar whereas the behavior of the high-loss modes has 
the opposite behavior with power of energy dissipation large and quality factor always small. 

5 Perturbation analysis of a system subjected to har- 
monic forces 

In this section we give an asymptotic description of the stored energy, power of dissipated 
energy, and quality factor for a harmonic solution v(t) = ve~ lwt of the system ([8]) in the 
high-loss regime > 1 subjected to a harmonic external force / (t) = fe~ luJt with nonzero 
amplitude / G H and frequency uj G IR. We will state our main results in this section but 
hold off on their proofs until Section El 

To begin we recall that according to fj46|) -f l48|) . assuming u is not in the resolvent set of 
the system operator A = Q — i/3B, there is a unique harmonic solution v(t) = ve~ luJt to 
the system (jSJ) with the harmonic force / (t) = fe~ luJt whose amplitude v is given by 

v = St(w)/ =i[ul- (fi - i/^H 1 /, (120) 
%(u)=i[uI-A(p)]-\ 

where 21 (u) is the admittance operator. 

As was introduced in Section [2~3| the stored energy U, power of dissipated energy W^s, 
and quality factor Q = Qf jU associated with the harmonic external force / (t) = fe~ luJt is 
given by the quantities 

U = \(v,v), W dis = (3(v,Bv), Q = \u,\JL = \u,\ jA~~~~r i (121) 
2 W dis (3 (v, Bv) 

where Q is said to be finite if W&i S ^ 0. 

For the results in the rest of this section we assume that /, u are independent of the loss 
parameter /3. 

The techniques of analysis in the high-loss regime ^> 1 differ significantly depending on 
whether the frequency u is an asymptotic resonance frequency or not. 

Definition 19 (nonresonance frequency) A real number u is called an asymptotic non- 
resonance frequency of the system (d)) provided u ^ pj, N B + 1 < j < N, otherwise it is an 
asymptotic resonance frequency. 

The usage of this terminology is justified by the following proposition: 

Proposition 20 Let u G IR. Then the admittance operator 21 (u) is analytic at (3 = oo if 
and only if u is an asymptotic nonresonance frequency. 

In this paper we will only consider the nonresonance frequencies. 
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5.1 Nonresonance frequencies 



In this section we state the results of our analysis of losses for external harmonic forces with 
asymptotic nonresonance frequencies in the high-loss regime. 

Recall from (TIT]) the Hilbert space H decomposes into the direct sum of orthogonal sub- 
spaces invariant with respect to the operator B > 0, namely, 

H = H © Hi 

where Hb = ran B, Hg = ker B are the loss and no-loss subspaces with orthogonal projec- 
tions Pb-i -Pg , respectively. It follows from this and the block representation of Q and B in 
(Fl8|) that £1 — A (j3), with respect to this decomposition, is the 2x2 block operator matrix 
in ( 152]) . namely, 



ojI - A ((3) = 
S 2 (u,f3) := ojI 2 - (fi 2 



5 2 (oj,/3) 
-6* 



-0 
Si (oj) 

Hi M : 



(122) 



cull — 



where Q\, Q 2 , and B 2 are self-adjoint operators, the latter of which has an inverse satisfying 
B^ 1 > 0. With respect to this block representation, the Schur complement of S 2 (oj, (3) in 
ojI — A((3) is the operator in ( I53p . namely, 



S 2 (u,P) = E 1 (oj)-Q*E 2 (u^r'Q, 



(123) 



whenever S 2 (oj, (3) is invertible. 

To simplify lengthy expressions we will often suppress the symbols oj, (3 appearing as 
arguments in the operators Si (oj), S 2 (oj, (3), S 2 (oj, (3). We now give the main results of this 
section for an asymptotic nonresonance frequency oj. 

Proposition 21 (admittance asymptotics) For (3 ^> 1, each of the operators Si (oj), 
E 2 (oj,(3), S 2 (oj,(3), and ul — A(f3) are invertible and the admittance operator 21 (oj) = 
i (ojI — A ((3))" 1 is given by the formula 



21 (oj) 







rs 2 > o " 




. 7 i 




s 2 l 





^2 



e*s^ 



^2 ©^2 ®*^2 "2 




h 



(124) 



Moreover, 21 (oj) is analytic at (3 = oo and has the asymptotic expansion 





21 (oj) 



o isr 1 



+ W^p- 1 + O (f3- 2 ) 



as (3 — >■ oo ; where 



1) = 


r B ? 




b^qi 


-1 


-1 






(sr 1 )*^^- 1 


( 




^esr 1 . 






h o " 




' B, 1 0" 






— 1 " 




_ (s r ye* h _ 









. J l 



(125) 



(126) 



In particular, this is a positive semidefinite operator, i.e., 

w { - 1] > 0. 
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Corollary 22 The operators 21 (u)* 21 (u), P B 2l (w)* 21 (w) P B; and 21 (u)* {3B% (u) are ana- 
lytic at (3 = oo. Moreover, they have the asymptotic expansions 



21* (w) 21 (w) 



+ 





o is^ecsrysr 1 



(127) 



-i 







and 




21* (u) /3P21 (u) = W^r 1 + O (/T 2 ) 



/T 1 + O (/T 2 ) , 
/T 2 + O (/T 3 ) , (128) 



:i29i 



as (3 — > oo. 



The following statements give our main results regarding the stored energy U, power of 
dissipated energy W^s, and quality factor Q = Qj^ associated a harmonic external force 
/ (t) = fe~ lult , where /, u are independent of (3, f ^ 0, and u an asymptotic nonresonance 
frequency of the system (jSJ). As we shall see the behaviour in the high-loss regime /3 ^> 
1 of these quantities is drastically different depending on whether the amplitude / has a 
component in the no-loss subspace H B or not, i.e., P B f ^ or Pgf = 0. 

Theorem 23 (quality factor) If P^f = then the stored energy U, power of dissipated 
energy Wdis, and quality factor Q are analytic at (3 = oo and have the asymptotic expansions 



b 2 2 + B^e (e 1 (u J y i y~ 1 (uj)- 1 e*^- 1 ] /) /r 2 + o (/r 3 ) 

W dls =(f,B 2 1 f) /3^ + 0(/3- 2 ), 



Q = M 



I (/, k" 2 + B?Q (S! (coy 1 )* S a (a;)" 1 0*5^ 



/ 



/r 1 + o (/r 2 ) 



(130) 
(131) 

(132) 



(f,B^f) 

as [3 —7- oo. Moreover, the leading order terms ofU and Wdis in these expansions are positive 
numbers, similarly for Q provided u ^ 0, and satisfy the inequalities 



- (7, [B 2 2 + B 2 l Q (E 1 (coy 1 )*-! (toy 1 G*B^\ / 



> 



sup 

1<3<N B 



{4}) "(/.^/)>i( i£ ^{o}) "(/,/)> 



Theorem 24 (quality factor) IfPgf ^ then the stored energy U and power of dissipated 
energy Wdis are analytic at (3 = oo and have the asymptotic expansions 



U = ~ (H! (u)- 1 P^f, S! (c)- 1 Pi/) + O (/T 1 ) 



(133) 



(/.W^r' + OC/T 2 ) (134) 

as /3 — >■ oo. In particular, the leading order terms in the expansions of U and Wdis are 
positive and nonnegative numbers, respectively. 
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Moreover, if Wdis ^ for /3 3> 1 then the quality factor Q has a pole at (3 = oo provided 
w^O. In particular, if f ^ ker W^ 1 ^ then it has the asymptotic expansion 

Q = M (/^(-i)/) ^ + C 1 ) ( 135 ) 

as /3 — > oo, whose leading order term is a positive number provided u ^ 0. 

Corollary 25 The operator W^~^ defined in U26\) has the N — Nb dimensional kernel 

ker W ( ~ l) = ker (P B + 0S X (c^ 1 F£) 

= {fi + f2eH:f 1 eH B - and f 2 = -BE, (a;)" 1 f,} . 

Corollary 26 (quality factor) The stored energy U , power of dissipated energy Wdi S , and 
quality factor Q have the following limits as (3 — )■ oo: 

lim U = ( ^ M - S. H- P B V) > i/ * 0, ( j 

i tfP£f = 0, 

lim = 0, (137) 

jfit— >oo 

/3->oo [ Pg j =0 oru = 0, 

where we assume in the statement regarding quality factor for the case Pg / 7^ £/ia£ Wdis ^ 
/or /3 ^> 1. In particular, if Pgf = or f ^ ker W^ -1 ) then Q is finite for (3 1 and we 
have the above limits for U, Wdis, and Q as (3 — > 00. 

6 Proof of results 

This section contains the proofs of the results of this paper. We split these proofs into 
two subsections. In the first subsection we prove the statements given in Section H] on the 
spectral perturbation analysis in the high-loss and low-loss regime for the system operator. 
In the second subsection we prove the statements in Section [5] on the perturbation analysis 
of losses for external harmonic forces. All assumptions, notation, and convention used here 
will adhere to that previously introduced in those two sections. 



6.1 Perturbation analysis of the system operator 

This purpose of this section is to prove the results given in Section HI We do so by considering 
in separate subsections the high-loss regime /3> 1 and the low-loss regime 0</9<l. 



6.1.1 The high-loss regime 

Proof of Theorem [5]. Let Q, 1 < j < Nb be an indexing of all the nonzero eigenvalues 
of B (counting multiplicities) where Nb = rank P. These eigenvalues are all positive real 
numbers since by assumption B > 0. 

We begin by extending the domain of the system operator by A {(3) = Q — i{3B, £ C 
Recall by our assumption fl is self-adjoint. Hence the operator (— i/3)" 1 A ((3) = P+(— i/3)~ Q 
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is analytic in /3 in a complex neighborhood of /3 = oo and if we use the substitution e = 
(—i/3) -1 then the operator eA (ie -1 ) = B + eVt is an analytic operator which is self-adjoint 
for real e. Thus by a theorem of Rellich |Bau85t p. 21, Theorem 1] we know for e G C with 
|e| <C 1, the operator B + eQ is diagonalizable and there exists a complete set of analytic 
eigenvalues Xj (e) and eigenvectors Xj (e) satisfying 

(S + eft) xj (e) = Xj (e) a;,- (e) , Xj (e) = Aj (e) , 1 < j < JV; (139) 
(xj (e) , x fc (e)) = 5j fc , for eel, 1 < j, k < N, 

where 5jk denotes the Kronecker delta symbol. In particular, the vectors Wj := Xj (0), 
j = 1, . . . , N form an orthonormal basis of eigenvectors for the operator B. Thus, after a 
possible reindexing of these analytic eigenpairs, we may assume without loss of generality 
that 

Bwj = Qwj, Xj (0) = Q > 0, 1 < j < N B] (140) 
Bwj = 0, Xj (0) = 0, N B + 1 < j < N. 

Denote the derivatives of these eigenvalues at e = by pj := X'j (0), 1 < j < N. Then it 
follows that they satisfy 

Pj = (wj, nwj) £l, 1 < j < N. 
Indeed, for e real with \e\ 1 by f |139[) we have 

= ( Xj (e) ,(B + en- Xj (e) I) Xj (e)) , l<j<N 

and the result follows immediately by taking the derivative on both sides and evaluating at 
e = 0. 

Now from these facts and recalling the substitution e = (— i/3) _1 that was made, we 
conclude that if /3 e C with > 1 then the operator A ((5) = (-i/3) (B + (-i/3)' 1 Q) is 
diagonalizable with a complete set of eigenvalues Q (/3) and eigenvectors lo, (/3) satisfying 

09) = (-i/3) Aj ((-i/3)" 1 ) , ^ 08) = Xj ((-i/3)" 1 ) , (141) 
A(p)w;(p)=C j (p)w j (fi), l<j<N. 

We will now show that the eigenpairs 

high-loss: Q i/3) , Wj 09) , 1 <j < N B ; 
low-loss: (/3) , ^' (/3) , iV B + 1 < j < ^, 

have the properties described in Theorem EJ 

We start with the high- loss eigenpairs. First, it follows from our results above that the 
eigenpairs Xj (e), Xj (e), 1 < j < N B are analytic for e G C with |e| < 1 and have the 
asymptotic expansions 

Xj (e) = Q + Pj e + O (e 2 ) , Q > 0, Pj G R, 
(e) =Wj + 0(e), 1 < j < 

as £ — >■ 0. This implies by fll4ip the high- loss eigenvalues have poles at (3 = oo whereas their 
eigenvectors are analytic at /3 = oo and they have the asymptotic expansions 

Q 09) = (-i/3) A, ((-i/3)- 1 ) = -ifiCj + p 3 +0 (/3" 1 ) , Q > 0, G K, 
Wj (/3) = Xj ((-i/3)" 1 ) = Wj + (/3" 1 ) , 1 < j < N B 
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as (3 — > oo, where the vectors Wj, 1 < j < Nb form an orthonormal basis of the loss subspace 
Hb = rani? and 

Bwj = Qwj, Pj = { w ji Qtij) , 1 < j < -^Vb- 

This proves statement (i) of this theorem. 

We now consider the low-loss eigenpairs. By the results above we know the eigenpairs 
\j (e), Xj (e), Nb + 1 < j < N are analytic for e G C with |e| <C 1 and have the asymptotic 
expansions 

X J (e)=p J e+ 1 -X';(0)e 2 + O{e 3 ), 
Xj (e) = Wj + x'j{0)e + 0{e), N B + l<j<N 

as e — > 0. Moreover, there is an explicit formula for the second derivative of these eigenvalues 

AJ(0) = 2 (4(0), Bx'j(0)) > 0, N B + 1 < j < N. 

Indeed, this follows from (11391) and (11401) since 

PjS + ^Aj(0)e 2 + O (e 3 ) = Xj (s) = (Xj (e) , Qxj (e)) £ + (^ (e) , fix, (e)) 
= (Wj, Qwj) e + (4(0), 5x;(0)) £ 2 + O ( £ ~ 3 ) 

for real e as e — >■ 0. Thus we can conclude from this and (1141R that the low-loss eigenpairs 
are analytic at = 00 and they have the asymptotic expansions 

(/?) = (-i/3) A, (H/3)" 1 ) = - id,/?- 1 + O (/T 2 ) , Pi G ffi, > 
^ = ((-i/?)" 1 ) = + + O (/T 2 ) , N B + l<j<N 

as /3 00, where the vectors u^, A?# + 1 < j < A r form an orthonormal basis of the no- loss 
subspace Hg = ker B and 

.Biyj = 0, pj = (wj, Qwj) , Wj' 1 ** = ix'j(0), 

- = -i§Aj(0) = -i (xJ(0),SxJ(0)) = -i (^i^) 

for A?# + 1 < J < A^. This completes the proof. ■ 

Proof of Corollary M Let Wj (/3), 1 < j < N B and Wj (/3), N B + 1 < j < N denote 
the high-loss and low-loss eigenvectors, respectively, given in the previous theorem. Then by 
our results we know that the zeroth order terms in their asymptotic expansions must satisfy 
PbWj = Wj, Pj^Wj = for 1 < j < N B and PbWj = 0, P^Wj = Wj for N B + 1 < j < N since 
H = H B ® and Pb, Pb are the orthogonal projections onto H B , H^, respectively. The 
proof of this corollary now follows. ■ 

Proof of Proposition Let (j ((3), 1 < j < N be the high-loss and low-loss eigenvalues 
given in Theorem |5j Then as described in its proof we can extend the domain of the 
system operator A (/3) = Q — i(3B, f3 G C and these eigenvalues can be extended uniquely to 
meromorphic functions in a neighborhood of /3 = 00 whose values are eigenvalues of A (/3). 
Moreover, using the same notation to denote their extensions, it follows from (I139P and (I14ip 
that these functions satisfy 

OW=Ci(-0)> i<j<N 
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for all p e C with \(3\ > 1. This implies § (Q (J3) + Q (~P)) and ^ (£,■ (P) - Cj (-/?)) are 
even and odd functions, respectively, and meromorphic at p = oo and equal Re Q and 
ImCj respectively, for real (3. In particular, this implies their Laurent series in have 
only even and odd powers, respectively. The rest of the proof of this proposition now follows 
immediately by considering the real and imaginary part of the asymptotic expansions of the 
high-loss and low-loss eigenvalues given in Theorem |5j ■ 

Proof of Proposition [8} 1. &z 2. Here we will prove just the first two statements of 
the proposition. We come back to the proof of the third and fourth statements after we have 
proved Proposition [TT1 

Let j G {1, . . . , N} and (3 ^> 1. Well, since Q , Wj ((3) is an eigenpair of A ((3) then it 
follows from flUTl) and the fact B > that ImQ < 0. From this and Proposition since 
ImQ ((3) either as a pole or is analytic at P = oo, then either ImQ {[3) = or ImQ ((3) < 
0. If ImCj (13) = then it follows by (JSTJ that Bwj ((3) = and hence Cj (P)wj (p) = 
A ((3) wj ((3) = flwj (f3) which implies that Q (f3) = pj. Furthermore, if 1 < j < N B then we 
know that ImQ (f3) = -Qf3 + O as f3 ->■ oo with Q > and so ImQ (/3) < for p > 1. 
This completes the proof of the first two statements. ■ 

Proof of Proposition II 1L We begin by proving the first statement of this proposition. 
Recall that B2 = PbBPb\ Hb : Hb — > Hb where Pb is the orthogonal projection onto Hb- 

By Theorem [5] it follows that Q, tbj, 1 < j < Nb are a complete set of eigenvalues and 
eigenvectors for B 2 with these eigenvectors forming an orthonormal basis for Hb- As these 
eigenvalues are all positive this implies B 2 > and, in particular, it is invertible and its 
inverse satisfies B^ 1 > 0. This completes the proof of the first statement. 

Next, we prove the second statement of this proposition. Well, by Theorem we know 
that the vectors Wj, Nb + 1 < j < N form an orthonormal basis for the no-loss subspace 
Hq. By definition Pg is the orthogonal projection onto implying PgWj = Wj, for 

j = Nb + 1, . . • , N. And hence, since f2i = PgQPM ± : — > which is obviously a 

B 

self-adjoint operator because Pg and Q are, we have QiWj = P^Qwj, j = Nb + 1, . . . , N. 
Therefore this, the fact P^P = 0, and Theorem [5] imply 

tttWj = P^VtWj = lim P^ (D, - ipB) Wj (p) = lim P^Q (p) Wj (p) = pjWj 

for j = Nb + 1, . . • , N. This completes the proof of the second statement. 

Finally, we prove the third and final statement of this proposition. Recall the operator 
6 = PeOPg-l x : -> H B - We now show that 

Qwj = iB 2 P B w { r l \ N B + 1 < J < N. (142) 
Well, by d92D, (USD, and <$T§ of Theorem [5] it follows that if N B + 1 < j < N then 

PjWj + O (p- 1 ) = Q (P) wj (P) = A (P) Wj (P) 
= (ft - ipB) Wj + /T 1 (ft - ipB) w { r l) + O (p- 1 ) = ttwj - \Bw { r 1] + O (p- 1 ) 
as P — > 00. Equating the zeroth order terms we conclude that 

PjWj = Qtij — iBwj ^. 
Applying Pg to both sides of this equation we find that 

= P B nwj - iP B Bw { r 1] = (Pb^Pb) ^j - i [PbBPb] P B w { r 1] = Q Wj - iP^P^^ 
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which proves the identity (j!42p . 

Therefore by (157)1 of Theorem [5l (I142p . the facts P B is the orthogonal projection onto 
H B , P B BP B = B, B 2 = P B BP B \ Hg , and since B 2 l is self-adjoint we conclude that 

-idj = -i (w^ x \Bw^ = - {p B w { f x) ,iB 2 P B w { r X) ) = - (p B w { r x \&w>) 

= - ^B 2 1 B 2 P BW y i \ew^ = -i (iB 2 P B w^ x) , B^Gwj) = -i (Gwj^^Gwj) 

= -i (w j ,&*B 2 1 &w J ) . 

This proves the final statement and hence the proof of the proposition is complete. ■ 

Proof of Proposition [§} 3. & 4. We now complete the proof of Proposition [H] by prov- 
ing the last two statements. Let j £ {1, . . . , N}. We begin by proving the third statement. 
Suppose that dj = 0. Then by Proposition [TTI we know that = — i (wj, Q*B 2 ~ 1 Qwj) and 
B 2 X > 0. This implies Qwj = 0. By Theorem [5] we have Wj £ H B so that P B Wj = Wj and 
hence 

= Gwj = P B nP^Wj = Pettwj. 
From this, the fact P B + P B = I, and by the first statement in Proposition [TTJ it follows that 

{Ivjj = P B VlWj = P B VtP B Wj = QiWj = PjWj. 

Thus we have shown if dj = then ibj £ ker (pjl — Q). We now prove the converse. Suppose 
Wj £ ker (pjl — Q). This hypothesis and the fact Wj £ imply 

Qvjj = P B QP B Wj = P B Vtujj = pjP B Wj = 0. 

Therefore by Proposition [TTI we conclude dj = ^Wj,Q*B 2 ~ 1 Qwj) = 0. This proves the third 
statement of Proposition [HI 

Finally, we will complete the proof of Proposition [HJ by proving the fourth statement. 
Suppose Wj (jL ker (O) but pj = and dj = 0. Then it would follow from the third statement 
of this proposition that Wj £ ker (Q), a contradiction. Thus we have shown if Wj ^ ker (Q) 
then pj 7^ or dj ^ 0. We now prove the converse. Suppose pj ^ or dj ^ 0. If 
Wj £ ker (Q) then 

pjWj = £l\Wj = P B VtWj = 0. 

This implies pj = and hence Wj £ ker (pjl — Q). By the previous statement in this 
proposition this implies dj = 0. This yields a contradiction of our hypothesis. Therefore 
we have shown if pj ^ or dj ^ then Wj £" ker (Q). This proves the fourth statement and 
hence completes the proof of this proposition. ■ 

Proof of Corollary 1121 The proof of this corollary is straightforward and follows from 
Theorem [5] and Proposition [HJ ■ 

Proof of Proposition [T3l For high-loss and low- loss eigenvectors Wj 1 < j < N we 
have (wj , Wj (/3)) = 1 +0 as (3 — > oo. From this it follows that the energy has the 

asymptotic expansions 

U [wj 09)] = \ (Wj , Wj (/?)) = i + O (r 1 ) , 1<J<N 

as (3 — > oo. Now by Proposition [7J we know the imaginary parts of high-loss and low-loss 
eigenvalues Q 1 < j < N have the asymptotic expansions 

lmQ(f3) = -( 3 f3 + 0(f3~ 1 ), Q > 0, 1 < j < N B ; 
Im Q (/3) = -dj/3- 1 + O (/T 3 ) , dj > 0, N B + 1 < j < N 
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as (3 — > oo. Thus by these asymptotic expansions we conclude the power of energy dissipation 
of these eigenvectors by the formula fl9T|) have the asymptotic expansions 

w^tov® = -2 mow utov® = { JfS + g^) if ^^'Ifl'jv 

as /3 — > oo. The proof of the proposition now follows from these asymptotics. ■ 

Proof of Proposition 1141 By the formula (l9Tj) for the quality factor of an eigenvector 
of the system operator and by Proposition [7J it follows for (3 3> 1 we have Q [wj (/?)] is finite 
if and only if Im(j ^ in which case 

By Proposition [7] we know that Re (j and Im Q (f3) are either analytic or have poles at (3 = 
oo whose Laurent series expansions contain only even and odd powers of respectively. 
In particular, this implies either KeQ {(3) = or the sign of Re(j {(3) does not change for 

1 and |ReCj = Re(j (/?), for /3 3> 1 if this sign is positive or if this sign is negative 
|ReCj 0#)| = — ReCj for (3 ^> 1. It follows from this and Proposition [7J that we have the 
asymptotic expansions 

|ReO 03) | = \p,\ + O (f3- 2 ) , ImQ ((3) = -Qf3 + O (/T 1 ) , 1 < j < N B ; 

|ReO(/?)| = \pj\ + o (/r 2 ) , imc, (£) = -^/r 1 + o (/r 3 ) , Afe + 1 < j < n 

as (3 — > oo. It also follows that if Im(^ ^ for /3 ^> 1 then Q ((3)} is the product of 
two meromorphic functions one of which is even and one of which is odd in a neighborhood of 
(3 = oo which implies it has a Laurent series expansion at (3 = oo containing only odd powers 
of and, in particular, lim^oo Q [wj (/?)] = or oo. These facts and the fact Q > for 

1 < j < Nb implies we have the asymptotic expansions 



Q [wj 



-l\ReQW)\ j ^(3- 1 + 0((3- 3 ) if 1 < j < N B , 



as (3 — > oo. The proof of this propositions now follows from this. ■ 



6.1.2 The low-loss regime 

Proof of Theorem 1151 Let Uj, 1 < j < N be an indexing of all the eigenvalues of Q 
(counting multiplicities). The proof of this theorem is similar in essence to the proof of 
Theorem [51 We begin by extending the domain of the system operator by A (/3) = Q — i(3B, 
(3 G C. Using the substitution e = —i(3 then the operator A (ie) = Q + eB is an analytic 
operator which is self-adjoint for real e. Thus by a theorem of Rellich |Bau851 p. 21, Theorem 
1] we know for e G C with |e| <C 1, the operator Q + eB is diagonalizable and there exists a 
complete set of analytic eigenvalues (e) and eigenvectors Xj (e) satisfying 

(fi + eB) Xj (e) = Xj (e) Xj (e) , X~ (I) = Xj (e) , 1 < j < (143) 
(e) , x k (e)) = 5 jk , for £GK, 1 < j, k < N, 
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where 5jk denotes the Kronecker delta symbol. In particular, the vectors Uj := Xj (0), 
j = 1, . . . , N form an orthonormal basis of eigenvectors for the operator Q. Thus, after a 
possible reindexing of these analytic eigenpairs, we may assume without loss of generality 
that 

Quj = cOjUj, Xj (0) = Uj, 1 < j < N. 

Denote the derivatives of these eigenvalues at e = by Oj := A^ (0), 1 < j < N. Then it 
follows that they satisfy 

Oj = (uj, Buj) el, l<j<N. 
Indeed, for e real with \e\ <C 1 by ( I143p we have 

= ( Xj (e) ,(Q + eB- Xj (e) I) Xj (e)) , 1 < j < N 

and the result follows immediately by taking the derivative on both sides and evaluating at 
e = 0. 

Now from these facts and recalling the substitution e = —i/3 that was made, we conclude 
that if ft G C with \ft\ 1 then the system operator A = Q + (—i/3) B is diagonalizable 
with a complete set of eigenvalues Q ((3) and eigenvectors Vj (f3) satisfying 

09) = A, (-i/3) , vj 09) = (-^) , (144) 
A{fi)v j {p) = Q{p)v j {p), l<j<N. 

We will now show that these eigenpairs have the properties described in Theorem [15J First, 
it follows from our results above that the eigenpairs Xj (e), Xj (e), 1 < j < N are analytic for 
e G C with |e| <C 1 and have the asymptotic expansions 

Xj (e) = Uj + OjE + O (e 2 ) , Uj G M, cr^ > 0, 
Xi ( e ) = Uj + O (s) , 1 < j < N 

as e — > 0. This implies by f)144p these eigenvalues and eigenvectors are analytic at ft = 
and they have the asymptotic expansions 

Q (ft) = Xj (-i(3) = Uj - Wjft + O (ft 2 ) , Uj G R, Uj > 0, 
^ (ft) = x, (-i/3) = Uj + O (ft) , l<j<N 

as ft — y 0, where the vectors Uj, 1 < j < Ng form an orthonormal basis of eigenvectors of Q 
and 

Quj = UjUj, Oj = (uj, VLuj) , 1 < j < N. 

This completes the proof. ■ 

Proof of Corollary 1161 This corollary follows immediately from the asymptotic expan- 
sions (11141) of Theorem [15] and the formulas from ( J91~l) . ■ 

Proof of Corollary I17L The proof of this corollary is similar in essence to the proof 
of Proposition [TU As in the proof of Theorem [15] we can extend the domain of the system 
operator A (ft) = Q — iftB, ft G C and the functions Q (ft), 1 < j < N can be extended 
uniquely to analytic functions in a neighborhood of ft — whose values are the eigenvalues 
of A (ft). Moreover, using the same notation to denote their extensions, it follows from (I143P 
and (11441) that these functions satisfy 

GW = G(-/3), l<J<iV 
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for all G C with \0\ < 1. This implies § (/?) + Cj (~P)) and ^ (£,■ (0) - Cj (-£)) are 
even and odd functions, respectively, are analytic at = 0, and equal ReC,- (0) and ItnQ (0), 
respectively, for real 0. In particular, this implies their Taylor series in have only even and 
odd powers, respectively. And this implies as in the proof of Proposition [TH that \ReQ (0)\ 
is analytic at = and its Taylor series in has only even powers. Now by the formula in 
(I9T1) we know that if Im(j (0) = for G R with -C 1 then the quality factor is given by 

Q[VAIj)l 2 Im0(/3) ' 

and hence can be extended to a function meromorphic at = whose Laurent series in 
contains only odd powers. Using these facts and the asymptotic expansions in flll4j> of 
Theorem IT51 we conclude 



Q[vj 



l|ReCi(/3)| _ 1 \uj\ + O(0 2 



2 lm(j(0) 2-a 3 + O(0 3 ) 
1 U V + 0(/3) 



2 a, 



as /3 — > oo, provided cr^ 7^ 0, in which case lmig^ Q [ v j (0)] = or 00 depending on whether 
uij — or uij ^ 0, respectively. This completes the proof. ■ 



6.2 Perturbation analysis of a system subjected to harmonic forces 

This purpose of this section is to prove the results given in Section HJ 

Proof of Proposition 1201 Let u 6 R. If u is an asymptotic nonresonance frequency of 
the system ([H]) then Proposition [2H which we have proved in this section below, tells us the 
admittance operator 21 (u) is analytic at = 00. So in order to complete the proof of this 
proposition we need only prove the admittance operator 21 (a;) cannot be analytic at = 00 
if uj is an asymptotic resonance frequency of the system ([8]). Thus suppose uj = pk for some 
k G {Nb + 1, . . . , iV} but 21 (uj) was analytic at — 00. Then it is continuous and so by 
flS2D, (JHSD this implies 

= lim [-i (ul - C fc (0)) 21 (u) w k (0)} = lim w k (0) = w k ^ 0, 

a contraction. This contradiction proves the proposition. ■ 

For the rest of the proofs in this section the symbol uj whenever it appears will mean an 
asymptotic nonresonance frequency of the system ([8]), i.e., uj G M. and uj 7^ pj, Nb + 1 < 
j < N. Also, without loss of generality we may extend the domain of the system operator 
A (0) = Q — i0B, G C. Similarly we extend the domains of the operators S 2 (uj,0) = 
ul 2 - (fla - i0B 2 ), G C and S 2 (uj, 0) = S x (uj) - B*E 2 (uj, 0)' 1 6, G C provided E 2 (uj, 0) 
is invertible. For the rest of the proofs in this section we use these extensions. Also whenever 
it is convenient we will suppress the dependency of these operators on the symbols uj, 0. 

Proof of Proposition 1211 We begin by proving for \0\ 3> 1, the operators Si (u), 
E 2 (uj,0), S 2 (u,0) are invertible with E 2 (u,0)~ 1 , S 2 (uj,0), and S 2 (u,0)~ 1 analytic at 
— 00. By Proposition [TT] it follows that the spectrum of the operator Q± is the set 
{pj : N B + 1 < j < N} and since uj ^ {pj : iV^ + 1 < j < N} then the operator Hi (uj) = 
ujI\ — Qi is invertible. Now recall the well-known fact from perturbation theory that if 
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T is an linear operator on H such that in the operator norm ||T|| < 1 then I — T is in- 
vertible and the series Y^^Lo^ converges absolutely and uniformly to (J — T)~ . This 
implies if T is an operator- valued function which is analytic at (3 = oo then it has an 
asymptotic expansion T {(3) = T + + O {(3~ 2 ) as (3 — > oo and if T is invertible then 

T ((3) = T (/ — Tq 1 [T {(3) — T ]) is invertible for 3> 1 and its inverse T is ana- 

lytic at /3 = oo with the asymptotic expansion T(/3) _1 = Tq 1 — T^TiT^ 1 (3~ l + O {(3~ 2 ) 
as /3 — > oo. We will use some of these facts now to prove S 2 (a;, (3), S 2 (u>, (3) are invertible 
for \(3\ ^> 1 and S 2 (w,/3) _1 , ^(w,/?), and S2 (u>, /3)~ 1 are analytic at /3 = 00. Well, the 
function T {(3) = i/3 _1 H 2 (u>, (3) is analytic at (3 = 00 and its limit as /3 — > 00 is an invertible 
operator since 

lim i/T 1 ^ (C, (3) = lim i/T 1 [C/ 2 - (fi 2 - i/?S 2 )] = £? 2 > 0. 

/?— >oo /3->oo 

This implies S 2 (w, /3) = — i/3T is invertible for ^> 1 with the inverse H 2 /3) 1 = 
i/3 _1 T(/3) _1 analytic at /3 = 00 having the asymptotics 



H 2 (w, = -\B 2 1 (3- 1 + BJ 1 (C/ 2 - fi 2 ) ^/T 2 + O (/T 3 ) 
as /3 — ?• 00. From which it follows that the operator 

S 2 {u, (3) = E x (w) - 0*H 2 (w, = 2! (u) (h - S x (a;)" 1 0* ~ 2 (w, 0) 



(145) 



is well-defined and invertible for \/3\ 3> 1 as well as it and its inverse, S 2 (u, fS)~ , are analytic 

at (3 = 00. 

Now we prove — A ((3) is invertible for \/3\ 3> 1 and 21 (w) = i (u;I — A {f3))~ 1 is analytic 
at (3 = 00. First, by the 2x2 block operator matrix representation of ul — A {(3) from 
(11221) and since E 2 (u,(3) is invertible for \(3\ ^> 1 then, as discussion in Appendix [7] on 
the Atiken block diagonalization formula fl 1 49 [) — (11511) . the operator admits for \(3\ 3> 1 the 
Frobenius-Schur factorization 



ul - A (/?) 



S 2 -0 " 




-0* Si 





W ^ 2 





h 



s 2 
5 2 



h 




/1 



Furthermore, this implies for f3 ^> 1 that since S 2 (u), /3) is invertible then a;/ — A (/3) is 
invertible and 



2t(w) 



" / 2 s 2 x e " 




[sj 1 " 




/ 2 " 


. J i 




sv 1 . 




. e*s 2 1 h _ 



(146) 



"2 + "2 ©* ,=; 2 "5 



«5 2 - 1 e*s 2 1 



which proves formula (11241) . From this formula and the fact both H 2 (u, (3) 1 and S 2 (u>, (3) 1 
are analytic at (3 = 00 we conclude that 21 (w) is analytic at f3 = 00. Moreover, this formula, 
(I151l) - (ll54p in Appendix [3, that fact(Si (oj) *) = Hi (uj) 1 , and (11451) imply we have the 
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asymptotic expansion 



21 M 



J 2 



— 1 ~ J —2 





o i-r 1 



"2 

-lcx*^— 1/. 



+ W { - X) P~ X + O (/T 2 ) , where 

2 



O (/T 2 ) 
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as j3 — 7- oo. As PT 1 > it follows immediately from this block operator representation that 
PJA -1 ) > 0. This completes the proof of the proposition. ■ 

Proof of Corollary 1221 As the the operator 21 (oj) = i (ojI — A is analytic at 

/3 = oo and since A (/?)* = A(-f3) for real then 2t(w)* = -i (ul - A for /3 

real which is clearly analytic at /3 = oo. This implies each of the operators 21 (a;)* 21 (oj), 
P b $1(u)*QI(u)P b , and 21 (u)* S21 (w) is analytic at £ = oo. The fact 21 (w)* /3P21 (w) is 
analytic at (5 = oo follows immediately from our proof below that linig^oo 21 (a;)* -B21 (a;) = 0. 
Now the asymptotic expansion of 21 (u)* 21 (u) as — > oo in fl 1 2T[) follows immediately from 
the asymptotic expansion of 21 (u) in (j!25p and the definition of the real part of an operator 
T as ReT = \{T + T*). The proofs of the asymptotic expansions of Pb21 (oj)* 21 (a;) Pb 



and 21 (a;)* /3P21 (a;) in (11281) and (11291) are similiar, using the asymptotic expansion (11251) for 
21 (a;) and the fact that P B and B = P B BP B , with respect to the direct sum H = H B © H^, 
are the block operators 

(148) 
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+ B 2 1 Q(~^y~^ 1 Q*B 2 1 




It follows from ffT26j) and ffT25|) that 

P B (W^)*W^P B -- 

(W^Y BW^ =W { ~ 1 \ 
P B 2l (oj) = PsW^p- 1 + O (/T 2 ) , 21 (oj) P b = W { - l) P B p- 1 + O (/T 2 ) , 
P B Ql (oj)* = P B + O (/T 2 ) , 21 (oj)* P b = (W^yPsp- 1 + O (/T 2 ) 

as (5 — > oo. Thus from these facts it follows that 

P B 2l (oj)* 21 (oj) P b = P B (W { - 1] y W { - l) P B f3~ 2 + O (/T 3 ) 

21 (oj)* P2t (to) =%(oj)* P B BP B % (oj) = (W^)* PbBPbW^P' 2 + O (/T 3 ) 

= w^r 2 + o (r 3 ) 

as (3 — > oo. Therefore (I128f) is proven and from the fact 21 (a;)* P21 (a;) is analytic at f3 = oo 
with the asymptotic for it just derived, it follow that 21 (oj)* /3P21 (oj) is analytic at j3 = oo 
and has the asymptotic expansion (I129p . This completes the proof. ■ 

Before we now proceed to prove Theorems [23] and [24] we will find it convenient to first 
prove Corollary [25] 
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Proof of Corollary 1251 From the block operator factorization of the operator ^ in 
( I126P with respect to the direct sum H = H B © H B it follows immediately that W^~^ f = 
if and only if P B f + GHj (u)' 1 P£f = 0, i.e., / G ker (P B + 0Si P^) . But from this 
direct sum and the fact Pb, P b are the orthogonal projections onto H B , H B , respectively, 
we can compute this kernel to conclude 

ker W { - 1] = ker (P B + QEy (co)' 1 P%) 
= {fi + f 2 eH:f 1 eHi and f 2 = -6~i (co)- 1 }\) . 

This representation of the kernel, the direct sum, and the fact that dim H B = N — N B implies 
the kernel has dimension iV — N B . This completes the proof. ■ 

Proof of Theorems [23] & 1241 Let / G H, f ^ 0. The stored energy U and power of 
dissipated energy W&s associated with the harmonic external force / (i) = fe~ lu)t by (j!20p 
and ( I12ip are given by the formulas 

U = i (/, 21 (u;)* 21 (w) /) , ^ dis = (/, 21 (w)* pm (cj) f) 

for > 1. In particular, by these formulas and Corollary [22] it follows that U and are 
analytic at = oo and by (11271) . (11291) have the asymptotic expansions 

^ dis = (z,^- 15 /) r x +o(p- 2 ) 

as /3 — )■ oo. In particular, by the fact W^~ l "> > 0, the leading order term for W<±is is a 
nonnegative number and if P B f ^ then the leading order term for U is a positive number. 

Now since W<n s is analytic at (3 = oo then either W^is = for (3 ^> or Wdis 7^ for 
/? 1. Hence by the definition in (11201) and (11211) of the quality factor Q = Qf iU it will be 
finite for /3 ^> 1 if and only if W^s ^ for /3 ^> 0, in which case it is given by the formula 
Q = \oj\ U/Wdis implying it is a meromorphic function at /3 — oo. For example, by the 
asymptotic expansions just derived it follows that if P B f ^ then it must have a pole and 
if P B f = then it must be analytic. In particular, if we have (f,W^~^f^ ^ 0, which is 
equivalent to / G kerW^ 1 ^, then it has the asymptotic expansion 

o - u U^^p^f^^p-f) 
- M (KW^T) + (1) 

as (3 — > oo whose leading order term is nonnegative and if P B f ^ then it is positive. 

Now we complete the proof of Theorem [23] Suppose P B f = 0. Then we have P B f = f. 
Then it follows from this and the block operator representation for Vt^ -1 ) and P B in (1 1 2 6 [) 
and (I148p . respectively, that PsW^'^ P B f = B^ 1 f ■ In particular, since B^ 1 > then 
(/, W^ 1 ^ /) = (f, P^V) > and hence from the statements in this proof above the quality 
factor Q is finite for /3 > 1 and is analytic at /3 = oo. From our discussion above and (I128P 
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we have the asymptotic expansions 

w dis = (/,s 2 -7) /3- x +o(r 2 ), 

17 = ~ (/, a (a;)* 21 (w) /) = ± (/, P B 2l (a;)* 21 (w) P B /) 

= \ (7, W + TVe (~ r y h^tV] /) /r 2 + o (/r 3 ) 

as /3 — > oo. Thus to complete the proof of Theorem [231 we need only prove the inequalities 

described in that theorem. First, since B^ 1 > then it has a positive square root B 2 2 > 0. 
Second, it follows from Theorem [5] and fl 103[) that for any u G Hb we have 



> min C~ l = f max I 

l<j<N B 3 \l<j<N B J J 



2 



Thus imply, with u = f, B 2 2 f, the inequalities 



i(/,B 2 -'/)>(maX B 4-) "(/,/), 
i(/,B-/) = ( i m- B 6) W'/) 



2 

These facts imply, since P^G (S^ 1 )* H^ 1 6*B 2 ~ 1 > 0, that 



i (7, [b 2 - 2 + s^e (~rT -r 1 ^- 1 ] 7) > i 0) 1 (7, V/) 



1 

2 \l<j<N B 



>-[max 0) (7,7) >0. 



This completes the proof of the theorems. ■ 

Proof of Corollary 1261 The proof of this corollary follows immediately from Theorem 
and Theorem! 



Appendix: Schur complement and the Aitken for- 
mula 



Let M 



M 



P Q 
R S 



(149) 
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be a square matrix represented in block form where P and S are square matrices with the 
former invertible, that is, H-P^ 1 ]] < oo. Then the following Aitken block- diagonalization 
formula holds Zhang Sec. 0.9, 1.1], [Aitkent p. 67 (4)], 



M 



' P Q~ 




I " 




' P 




' I P- l Q ' 


R S 




RP- 1 I 




S P 




/ 



(150) 



i.e., the Frobenius-Schur factorization of the block matrix M |Tret08t p. xiv], where the 
matrix 

S P = M/P = S- RP- l Q (151) 
is known as the Schur complement of P in M. The Aitken formula (j!50p readily implies 



M" 1 



P- 



I -P~ l Q 
I 

P- 1 + P- l QSp l RP- 1 







Sp 1 



(152) 



I 

-RP- 1 I 

-P-'QSp 1 

Dp 

In particular, for ||-P _1 || <C 1 and under the assumption of invertibility of the matrix S, 
formulas f II 5 1 j) and f 1 1 5 2 j) imply 



-Sp^RP- 1 



Sp 1 =[S- RP^Q] 1 = S- 1 [I - RP^QS- 1 ] 1 
= S- 1 + S^RP^QS' 1 + O {\\P~ 2 \\) , 



(153) 



M" 1 





o s- 1 



M' 1 = [M -1 ] x + O (||P~ 2 ||) , where 



(154) 



[M~\ 



P- 1 -P^QS- 1 
S^RP- 1 S- 1 + S^RP^QS- 



p- 1 

S- 1 + S~ l RP- l QS- 1 



+ 



p- 1 
s- 1 



-Q 
-R 



p- 1 
o s- 1 



8 Appendix: Quality factor for eigenmodes 

In this appendix we derive a simple and descriptive formula for the energy Z7[iu], power 
of energy dissipation [w], and quality factor Q [w], for any eigenmode w of the system 
operator A (/?) with eigenvalue ( which are the quantities 

U[w\ = ~(w,w) (155) 
W dis [w] = {w, (3Bw) (156) 
„ r , _ energy stored in the system . 

Q M = 27[ t— \ 157 ) 

energy lost per cycle 
= |ReC|#^L = |ReC| 



W dis [w] 1 sl f3(w,BwY 
where Re£ denotes the real part of the eigenvalue ( and Q [w] is finite if Wdis M ^ 0. 
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Lemma 27 If w is an eigenvector of the system operator A (ft) = Q — \[5B with eigenvalue 
( then 

(w,ttw) (w,0Bw) 

ReC = ^ r , ImC = -, r-. 158 

(w,w) (w,w) 

In particular, if j3 > then lm( <0. 

Proof. By assumption B > and by hypothesis we have A(0)w = (w, w ^ so that 
(w, w) 7^ and 

Re(w,A(P)w) _ lm(w,A(P)w) _ 

(w,w) ' (w,w) 

On the other hand, since 

\ [A {(5) + A (/?)*] = n, i L4 - A (/?)*] = 



this implies 



Re (w, A w) (w,Qw) 
ReC = - 



lw,w) (w,w 



= Im (w, A (g) w) = (w,pBw) < Q 
(w,w) (w,w) ~ 

This completes the proof. ■ 

Proposition 28 (quality factor) If w is an eigenvector of the system operator A (j3) = 
Q — i(5B with eigenvalue ( then the energy U [w], power of energy dissipation W^s [w], and 
quality factor Q [w] satisfy 

Re(U[w] = - (w,Slw), W dis [w] = (w, (3Bw) = -2 lm(U [w] , (159) 

1 |( W ,^)| = _l^eC| 
WH 2(w,(3Bw) 2ImC' 1 J 

and Q [w] is finite if and only if Im( ^ 0. In particular, the quality factor Q [w] = — 

is independent of the choice of eigenvector w since it depends only on its corresponding 

eigenvalue (. 

Proof. By the lemma it follows that W<n s [w] = (w, (3Bw) = — (w, w) Im ( = —2 Im (U [w] 
and Re = | (w,w)Re( = \ (w,Qw). This implies that Q [w] is finite if and only if 
lm( 7^ 0, in which case 

or 1 ,p n U[w] 1 \(w,Qw)\ 1 \Regw,w)\ l\Re(\ 
^ m 1 ^'WdisH 2(w,pBw) 2-lm((w,w) 2 lm( ' 

This completes the proof. ■ 
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